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Abstract 

This review presents the effective temperature notion as defined from 
the deviations from the equihbrium fluctuation-dissipation theorem in out 
of equihbrium systems with slow dynamics. The thermodynamic meaning 
of this quantity is discussed in detail. Analytic, numeric and experimental 
measurements are surveyed. Open issues are mentioned. 
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1 Introduction 



One of the core ideas of statistical mechanics is that equilibrium states can 
be accurately described in terms of a small number of thermodynamic variables, 
such as temperature and pressure. At present there is no equivalent framework 
for generic out of equilibrium macroscopic systems, and one is forced to solve 
their dynamics on a case-by-case basis. 

Out of equilibrium macroscopic systems are of many different kinds. An 
interesting class is the one in which the relaxation is slow - with observables 
decaying, say, as power laws instead of exponentially. Typical instances are 
coarsening phenomena and generic glasses, realized as molecular, polymeric or 
magnetic materials, among others. Another intriguing group is the one of non- 
equilibrium steady states in the weak drive limit. Examples are gently vibrated 
granular matter and weakly sheared super-cooled liquids. 

The quest for an approximate thermodynamic description of such systems or, 
to start with, the identification of effective parameters acting as the equilibrium 
one, has a long history that we shall not review in detail here. In contrast, 
we shall focus on the development of the effective temperature notion, that has 
proven to be a successful concept at least within certain limits that we shall 
discuss. 

About 20 years ago, in the context of weak-turbulence, Hohenberg and 
Shraiman [1] proposed to define an effective temperature through the departure 
from the fluctuation-dissipation theorem (FDT). However, neither a detailed 
analysis of this quantity nor the conditions under which such a notion could 
have a thermodynamic meaning were given in this reference. Later, Edwards 
mentioned the same possibility in the context of granular matter [2]. Again, 
the proposal did not catch the attention of the community. In none of the en- 
suing studies the importance of distinguishing different dynamic regimes was 
sufficiently stressed and the few numerical checks performed at the time gave, 
consequently, confusing results. More recently, similar ideas appeared in the 
glassy literature [3, 4]. In this field, the possibility of solving exactly a number 
of schematic models [5, 6] put the definition of the effective temperature, Teff, 
on a much firmer ground. The solution to these models' dynamics established 
the importance of reaching an asymptotic limit with slow dynamics and of at- 
tributing a value of the effective temperature to each distinct dynamic regime. 
The solutions demonstrate the relation with the phase space volume visited dy- 
namically [7] , they illustrate in many forms the relevance of Teff to heat transfer 
and equilibration, and they set some limits to the extent of applicability of the 
concept. These results opened the way to studies in a wealth of more realistic 
cases. Beyond pure phenomenological descriptions of observed phenomena used 
in the past that arc of limited predictive power, see e.g. [8, 9, 10], the Teff def- 
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inition based on fluctuation-dissipation relations is not ambiguous and allows 
for direct measurements. 

In this review we sift through the effective temperature notion trying to 
transmit the full allure of the idea. The literature on fluctuation-dissipation 
violations is immense and several reviews have already been published [11, 12, 
13, 14]. We refer to these reports when appropriate. The rest of the review is 
structured in five Sections. The next one introduces the definition of a number of 
observables and the effective temperature. It follows a discussion of the insight 
gained from the solution to mean-field glassy models in Sect. 3. Section 4 
presents the interpretation of the effective temperature definition as a bona 
fide temperature. Section 5 is devoted to a (non-exhaustive) description of 
numeric and experimental measurements performed so far with special emphasis 
on recent studies. Finally, in Sect. 6 we present some brief conclusions. 

2 Definitions 

2.1 Canonical setting 

In this manuscript we focus on the dynamics of a classical or quantum sys- 
tem coupled to a classical or quantum environment, typically in equilibrium 
at temperature T. In this sense the setting is canonical. The systems relax 
by transferring energy to the environment. Atypical dissipation, as realized 
in granular matter, is referred to in Sect. 4.8.2. Reference to some studies in 
isolated systems (microcanonical setting) is made in the quantum context, see 
Sect. 5.13. All along this review ks = 1. 

2.2 Time and length scales 

In equilibrium studies one rarely specifics the initial time or implicitly takes it 
to —CO. Out of equilibrium a reference moment has to be defined. Time zero 
is usually taken to be the instant when the sample is prepared in some special 
conditions. In these notes we focus on infinitely rapid quenches and we define 
t = as the moment when the working temperature is instantaneously reached. 

We deal mainly with macroscopic systems in which the number of degrees 
of freedom is assumed to be the most diverging quantity, N oo. In most 
of the theoretical and numerical calculations discussed in this text times are 
finite with respect to N. Some reference to times that scale with the size of the 
system will be made. 

2.3 Correlations and responses 

The out of equilibrium dynamics are explored by measuring observables that 
depend on one or several times after the preparation of the sample. Observables 
could be of two types: those describing free evolution and those associated to 
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responses to external perturbations. A multi-time correlation is the average 
over histories or initial conditions of several observables Oj , that are functions 
of the system's degrees of freedom, evaluated at different times: 

Co^02...oAtl,t2,■■■,tn) = (01(^1)02(^2) ••.0„(i«)) . (1) 

Oj(t) is any observable with, without loss of generality, zero mean (otherwise 
one simply substitutes Oj Oj — (Oj)). 

In systems with an energy function the impulse response functions are the 
averaged reaction of an observable to a perturbation that modifies the potential 
energy as V ^ V — hOn at a given instant t„. One such response, at linear 
order in the perturbation, reads 

(2) 

Others, in which the perturbation(s) is (are) applied at intermediate time(s) 
can also be defined [15]. In dynamical systems in which the time evolution of 
the degrees of freedom is not related to the gradient of an energy function, one 
can similarly compute the response to an instantaneous kick that translates the 
degrees of freedom at a chosen instant as x(i„) x{tn) + Sx{tn) [16]: 

R f, , , ^ (5(0l(tl)...0„_l(t„„l)) 
^Oi02...0„ (,Cl, 12, • • ■ , trij = SxiT^ ' 

Causality implies that both _R's vanish when t„ > tj with j = l,...,n — 1. 

In cases with quantum fiuctuations ones needs to distinguish the order in 
which operators appear within the averages. The symmetric and anti-symmetric 
two-time correlations are 

C[OuO,]Jtl,t2) = (Oi(ii)d2(t2)) ± (02(i2)Oi(ii)) (4) 

with {Oi{ti)Oj{tj)) = Ti-[Oi{ti)Oj{tj)po]/Trpo, po the initial density operator, 
and Oi{ti) in the Heisenberg representation. The linear response to H ^ H ~ 
h02 at t2 is 

= C[o^,o;]Ah,t2) eih~t2) (5) 

h=Q 

and the last identity is the Kubo formula valid in and out of equilibrium at 
linear order in h. 

The search for a link between dissipation and fluctuations started with Ein- 
stein's derivation of a relation between the mobility and diffusion coefficient of 
a Brownian particle. The former is induced by an external force, and it is a re- 
sponse, while the latter is due to the spontaneous mean-square displacement of 
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P ... , . J(Oi(ti)...0.-i(<„_i)) 

dh{tn) 



Ro^O.{tl,t2) ~ 5/,(i2) 



the particle's position, thus a correlation. Onsager's regression hypothesis, and 
Kubo's linear response theory, elaborated upon the idea of an existing relation 
between the two types of fluctuations near equilibrium. In the rest of this Sec- 
tion we present some relations between induced and spontaneous fluctuations 
that hold in and out of equilibrium. 

2.4 Fluctuation-dissipation relations 

Several fluctuation-dissipation relations (FDRs) between responses and correla- 
tions computed over unperturbed trajectories have been derived in recent years. 
We mention just two of these FDRs that hold in and out of equilibrium. 

In generic Langevin dynamics with white [17] or coloured [18] equilibrium 
noise, and Markov processes described by a master equation with detailed bal- 
ance [13], the linear response is related to a correlation function of the unper- 
turbed system which, however, is more complex than the overlap between the 
two observables involved in the response. A simple example is the Id white- noise 
Langevin equation 

mi(t) + ji{t) = F[x{t)] + , (6) 
that using the fact that 2"fTR^{t,t') = {x{t)^{t')) implies [17] 

Ti?.(M') = l^^^ - l^^^ - l{-{mxit')]-x{t')F[xm . (7) 

F[x] is any deterministic force - not necessarily conservative. The generaliza- 
tion to other observables or field-theories is straightforward. The extension to 
the non-Markovian case is more easily performed in the path integral formal- 
ism, by exploiting the transformation of the generating functional under time- 
reversal [18]. This formulation allows one to prove that these relation holds for 
Hamiltonian dynamics as well. 

In the case of an Ising variable, s = ±, the equivalent relation reads [13] 

TR,{t,t') = \^^%^ + \{s{t)Y.W) - s"]wW) ^ s"]) , (8) 

s" 

with w the transition probability with no perturbation applied. This expression 
can be recast in the same form as eq. (7) since (s) = (Y^g„{s — s")w{s — >■ s")) 
and the second factor in the last term plays the role of a deterministic force. 
Extensions to multi- valued discrete variables are simple. This type of FDR can 
also be generalized to intrinsically out of equilibrium systems with transition 
rates that do not obey detailed balance. 

Many relations of this kind have been derived in the literature. A number of 
authors re-expressed them as a sum of two contributions: the ones that reduce 
to the FDT under equilibrium conditions and the anomalous ones that allow 
for many intriguing interpretations (dissipated energy flux, etc.). We refer the 
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reader to a recent review and a few articles that summarize these ideas [19] that, 
ahhough certainly very interesting, are not the main focus of this article. 

An FDR for dynamical systems characterized by a state variable x, that 
reaches an invariant stationary measure, say p{x), was proposed by Vulpiani et 
al. [16] Using the response to an impulse perturbation of type (3) one finds 

R^it,t') = ~m^-^^^^)e{t-t'). (9) 

Once again the generalization to a multi-variable system or a field theory and 
to other observables is immediate. Another relation of this kind was recently 
proposed by several authors [20]: 

R^it,t') = -{.(t)^±£^p^),^, e{t - 1') . (10) 

The average is computed at zero field. In order to go beyond this formal expres- 
sion one needs to know the measure ph- Different assumptions {e.g., Gaussian 
form [16], uniform [21]) lead to different relations. 

All these FDRs become the fluctuation-dissipation theorem (FDT) under 
equilibrium conditions, as explained in Sect. 2.5. 

2.5 Linear fluctuation-dissipation theorem 

The FDT states that the decay of spontaneous fluctuations cannot be distin- 
guished from the one of forced fluctuations. Its classical and quantum formula- 
tion are explained in the following paragraphs. 

2.5.1 Classical 

The classical linear and self fluctuation-dissipation theorem (FDT) expresses the 
equilibrium thermal fluctuations of an observable O in terms of the available 
thermal energy, T, and the linear response of O to a vanishingly small applied 
field linearly coupled to itself. FDT is used, for example, to infer mechanical 
properties of soft matter from the fiuctuations in light scattering. Importantly 
enough, there is no condition on the scale at which the observables are defined, so 
it could range from the macroscopic to the microscopic. Indeed, if one monitors 
the dynamics of a classical system, either Newtonian or stochastic with detailed 
balance, and assumes that equilibrium with the bath has been reached, it follows 
easily that 

Roit,t')^^^,Co{t,t')e{t-t'). (11) 

(This equation can be inferred from eqs. (8) and (9). In equilibrium the last 
term in the right-hand-side of eq. (7) vanishes due to reciprocity and the two 
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first ones are equal. If p takes the canonical form eq. (11) follows immediately 
from eq. (9). An elegant proof for Langevin processes with generic multiplicative 
and coloured Gaussian noise is detailed in [18].) The integral of the response 
function over a time-interval running from to i, 

Xo{t,t^)^ dt'Ro{t,t') , (12) 

is a dc susceptibility that is easier to measure experimentally. The FDT implies 
a linear relation between xo and Co 

Xo{t,t^) = ^[Co{t,t) - Co{t,t^)] . (13) 

A parametric plot of Xo/Co{t,t) against Co/Co{t,t), constructed at fixed t^, 
and for increasing t — t^,, is a straight line with slope — 1/T joining (1,0) and 
(OjXeq = 1/2^)- The same result is obtained by keeping t fixed and letting 
vary from to t. A departure from the straight line (13) for any observable O 
signals a divergence from equilibrium. 

A straightforward generalization is to the case in which the monitored ob- 
servable, A, is not the same that couples to the perturbation, B: 

RABitX) = ^-^CabM eit~t') . (14) 

In the body of this review and without loss of generality we focus on the self 
correlation and response, we take C'oit,t) = 1, and we erase the sub- index O 
until discussing observable dependencies in Sect. 5. Equation (13) is then more 
compactly written as x(C) = (1 — C)/T. 



2.5.2 Quantum 

When quantum fluctuations are important one needs to take into account the 
statistics of the observables at hand. Let and 0^ be (bosonic or fermionic) 
annihilation and creation operators, respectively. In the Schwinger-Keldysh 
closed time-contour formalism apt to deal with real-time Green functions these 
are defined as 

ihG''\t,t') = irm'it')) = Tr [Tc cl,H{t,a) 0l,(i',6) eH(0,±)] , (15) 

a, & = ±. (j) is either the complex conjugate (for bosons) or the Grassmann 
conjugate (for fermions) of (j). (j)ii{t,a) denotes the Heiscnberg representation of 
the operator at time t on the a-branch of the Keldysh contour. gniO, ±) = q{0) 
is the initial density matrix (normalized to be of unit trace) and its location on 
the -|- or — branch is not important thanks to the cyclic property of the trace. 
In the grand-canonical ensemble an equilibrium initial density operator reads 
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g(0) oc exp {—/3(H — ^N)), where A'^ is the number operator commuting with H 
(in non-relativistic quantum mechanics) and /i is the chemical potential fixing 
the average number of particles. Tc is the time-ordering operator acting with 
respect to the relative position of {t, a) and {t', b) on the Keldysh contour [22, 23]. 
A linear transformation of the fields allows one to define Keldysh, retarded and 
advanced Green functions: 

G^^(t,t') = |[G+-(t,t') + G-+(t,t')], 
G^it,t') = [G+-{t,t') - G-+(t,t')] ©(i 
G^(t,t') = [G+-(t,t') - G-+{t,t')] Q{t' 

respectively. 

The quantum fluctuation-dissipation theorem reads 

-C 

ImG^iuj) (17) 

where Fourier transforms with respect to t — t' have been taken and C = ± 
in the bosonic and fermionic case, respectively. One can formally recover the 
canonical ensemble result by setting fi = 0. In particular, associating G^ to R 
and G^ to G in the notation of the previous classical paragraphs and setting 
/X = 0, in the bosonic case one has 

Rit ^t')='- r — e— (*-*') tanh ( ^\ C{uj) Bit - t') , (18) 

that in the pfi ^ Q limit reduces to the classical expression (14). 



-t'), (16) 



tanh ( ■^(^■^ ^ m) 



2.6 Eulerian and Lagrangian formalisms 

A class of over-damped Markov diffusion processes take an equilibrium form - 
with detailed balance - in the Lagrangian (co-moving) frame of the mean local 
velocity. The resulting stochastic process does not contain information about 
the non- vanishing probability current of the original Eulerian (laboratory) frame 
dynamics. In practice, though, it is difficult to compute the local velocity and 
the passage to the co-moving frame is hard to implement. Still, this observation 
allows one to derive a modified FDT in the Eulerian frame that includes the 
effect of the mean local velocity v as an additional term [24] : 

TRAB{t,t') = e{t - t')\dt.GAB[t,t') ^ {[v{t') ■^B{t')]A{t)))\ . (19) 

In equilibrium the second term vanishes and the FDT is recovered. 
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2.7 Fluctuation dissipation ratio 



Consider the following adiniensional object 

where we assumed t > t' . When equilibrium is reached Teff ~ T. Out of 
equilibrium this ratio can be used to define a, possibly two time dependent, 
effective temperature although the thermodynamic meaning of this quantity is 
not ensured a priori. Note the asymmetry between t and t' in the denominator 
when the dynamics is not stationary. Alternatively, experiments are usually 
done in the frequency domain in which w) = dt" e^^'^ R{t, t — t) and 
C{t,uj) = Jq dr e*'^'^ C{t,t — r) and a time- frequency dependent effective tem- 
perature is defined as 

The time-domain and frequency-domain definitions of Teff are not necessarily 
equivalent. If this quantity is to have a physical meaning they must coincide 
under reasonable conditions. We shall come back to this point in Sect. 4.3. 
The asymptotic values 

X°° = lim \im X{t,t') , Te^f=T/X°° (22) 

(or lim(^_>.o limt_j.oo -'^(i, w) in the frequency-domain) turns out to be useful in 
the description of critical quenches, see Sect. 5.3. 

It is clear that X can be computed under any circumstances but its interpre- 
tation as leading to an effective temperature will not always be straightforward. 
In the rest of this section we shall still use the name Teff; later in the review 
we shall discuss when this is justified and when it is not. 



2.8 Beyond linear-response 

Building upon the study of no n- linear responses in Langevin processes [15], 
Lippiello et al. presented a unified derivation of out of equilibrium fluctuation- 
dissipation relations (FDRs) to arbitrary order, for discrete (Ising or g-state) and 
continuous variables ruled by stochastic Markovian dynamics of quite generic 
type - single flip, Kawasaki, Langevin, etc. [25] These relations are conceptu- 
ally simple although pretty lengthy to write down. They are derived from a 
rewriting of the multiple variations of the evolution operator with respect to 
the perturbation strength. In short, one ends up with a sum of an increasing 
number of terms that are time-derivatives of higher order correlations. 

The generic FDRs can also be derived from a fluctuation principle that 
simply uses the detailed balance property of the transition probabilities [15, 
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18, 25]. The proof goes as follows. One first relates the probability of a path, 
{(/)} = ^(i), under the effect of a perturbation h{t) and conditioned to the initial 
value 4>{ti) = to the probability of the reversed path, {0^} = (j)R{t) = 4>{tii), 
under the reversed perturbation and conditioned to 4)R{ti) = 4>{tF) — (f'F- 

e-'^^^'' ''^'^'^'^=P[{M;^F,{hR}] . (23) 

From the average over all paths weighted with a generic distribution of initial 
conditions, P{(f>i), the expansion in powers of h{t) yields all non-linear FDRs. 
This route can be followed for discrete and continuous variables as well. 

These proofs do not use any equilibrium assumption and the FDRs hold 
in full generality (for aging systems, in non-equilibrium steady states,...). In 
equilibrium the use of the Onsager relation and the stationary property allows 
one to simplify the FDRs considerably. Still, as soon as the linear regime is 
left, the non-linear impulse responses are linked to sums of several derivatives 
of correlations. Interestingly enough, beyond linear order the FDRs depend on 
the microscopic dynamic rule not only out but also in equilibrium. 

The interest in the generalized non-linear FDRs is at least three-fold. First, 
they allow one to develop efficient algorithms to compute linear and non-linear 
responses without applying any perturbation. Second, they can be used to 
search for growing dynamic correlation lengths in glassy systems, a field of active 
research [26] . Third, and more importantly for the purposes of this review, they 
provide a way to further test the consistency of the effective temperature notion 
that, in the non-linear relations between responses and correlations, should play 
essentially the same role as in the linear ones. 

To our knowledge, quantum FDRs of this kind have not been derived yet. 

3 Insight 

Fully connected disordered spin models provide a mean-field description of 
glassy phenomenology. They capture many features of glassy thermodynamics 
and dynamics, at least in an approximate way. They have static and dynamic 
transitions at finite temperatures, Ts and Td, that do not necessarily coincide. 
In a family of such models aimed to describe fragile glasses with the random first 
order transition scenario (RFOT) [27, 28, 29] T, corresponds to an entropy crisis 
realizing the Kauzmann paradox. Below an infinite system does not come to 
equilibrium with the environment [5, 6] and relaxes out of equilibrium. Replica 
tricks [30] and the Thouless- Anderson-Palmer approach [31] give us a handle to 
understand equilibrium and metastable states as well as thermodynamic prop- 
erties. A dynamic treatment completes the picture explaining their relaxation 
and the relation with the free-energy landscape. Importantly enough, a quasi- 
complete analytic solution exists and a clear definition of many concepts that 
are not as sharply grasped in finite-dimensional cases - such as metastable states 
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- can be given. We do not intend to give here a full account of the behaviour of 
these models; several reviews exist [32, 33] already. We simply highlight some 
of their properties and their relevance to the effective temperature notion. 

3.1 Slow relaxation and convergence 

The energy density as well as any other intensive observable that depends on 
just one time approach a finite value: 

-oo < Ooo = lim lim (0(t)) < +oo . (24) 

With this order of limits, the intensive quantity Ooo is not necessarily equal to 
the equilibrium value Oeq- Below Trf the approach is achieved slowly, typically 
with a power law. If time is let scale with N a different dynamic mechanism 
sets in and equilibrium is eventually reached, limAr^.oo limt^oo(0(t)) = Oeq- 

3.2 Separation of time-scales 

Below their dynamic critical temperature, T^, mean-field glassy models relax in 
several time-scales (see [6] for a precise definition). The correlation and linear 
response admit a decomposition: 

C{t,t') = C^°\t,t') + C^^\t,t') + C^''\t,t') + ... (25) 
R{t,t') ^ R^'^^t,t') + R'^'\t,t') + R^^\t,t') + ... (26) 

The first of these scales, usually labelled ST, remains finite as i' oo, while the 
others diverge with t' . In mean-field models these time scales become infinitely 
separated as i' — > cxd due to their different functional dependencies on t' . This 
means that when one of the terms varies the other ones are either constant 
or have already decayed to zero. In particular, in the ST scale the correlation 
function decays from its equal-times constant value to the Edwards-Anderson 
parameter q^A that depends on bath temperature and all coupling constants. In 
this scale the system follows the rules dictated by the bath and the dynamics are 
stationary, while in the other scales interactions are relevant and the relaxation 
is more complex. Each term contributing to the correlation scales as 

and h^°\t) = e*/^*"' = /i'^(t) with r^o) finite and f^\oo) = q^A- Loosely, one 
can define a time-scale r'*' from 

- 1 + 7(7)7^ with r«(0 = [dlnh^Ht')/dt'r ■ (28) 
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One of the main ingredients in the solution to mean-field glassy models in 
relaxation [5, 6] is the fact that eq. (13) does not hold below Td- This does 
not come as a surprise since the equilibrium condition under which the FDT is 
proven does not apply. What really comes as a surprise is that the modification 
of the spontaneous/induced fluctuations relation takes a rather simple form: 

i?W(i,0 = _l^|.cW(t,t') 0(^-0 z = 0,l,... (29) 

in the t' ^ oo limit. The temperature of the external bath is the one of the 
stationary regime and the conventional FDT holds, T^^ = T^^^ = T. We discuss 
the reason for this in Sect. 4.8. In all other regimes the T(')s are different from 
each other and from T. If one keeps t' finite the crossover from one scale to 
another is smooth and a generic relation like (20) applies. 

So far we highlighted three characteristics: the slow relaxation, the fact that 
the energy density is bounded from below, and the separation of time-scales. 
They are all fundamental for finding an effective temperature with a thermal 
significance. We shall come back to these properties in Sect. 5. 



3.3 The parametric construction 

A very instructive way to study the deviations from FDT is to construct a 
parametric plot of the integrated linear response against the correlation at fixed 
tu,, varying t between ty^, and infinity. In the long tu, limit such a function 
converges to a limiting curve: 

lim x{t,t^) = / dC , (30) 

where the effective temperature [3], Teff(C'), is a function of the correlation C. 
The main aims of this article are to justify the name of this function and to 
give a state of the art review of its measurements in model systems, numerical 
studies and experiments (see [11] for a previous survey of measurements of 
FDT violations with numerical methods and [14] for a recent discussion on 
effective temperatures from a different perspective). The convenience of this 
parametrization is that it makes it possible to compare to replica calculations 
and, on a more physical level, it allows for a better interpretation of fluctuations. 

Let us summarize the behavior of the asymptotic x(C') function in different 
mean-field glassy models relaxing slowly with bounded dynamics after a quench 
from the disordered into the ordered phases. In all cases there is a regime 
t — tw = 0(1) in which x{(^) is a straight line of slope — 1/T, joining (1,0) 
and {qEA{T), [1 — qEA{T)]/T). A free-energy interpretation of this regime will be 
given in Sect. 4.8. At this point the straight line breaks and the subsequent 
behaviour depends on the model. Three families have been identified. 
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Models describing domain growth such as, for example, the ferromagnetic 
0{N) model in d dimensions in the large N limit. There is only one additional 
time-scale in which C decays from (Jea to 0. The asymptotic parametric plot for 
C < qEA{T) is flat. The susceptibility x gets stuck at the value [1 — gEA(T)]/T 
while the correlation C continues to decrease towards zero. The same result 
holds for the Ohta-Jasnow-Kawasaki approximation to the Xcj)'^ model of phase 
separation. In agreement with the fact that temperature is basically an irrele- 
vant parameter in (at least clean) coarsening the effective temperature diverges 
in the full low-temperature phase, T^^^ — >■ oo. 

Models of random- first- order-transition type (RFOT), that provide a mean- 
field-like description of fragile glasses. Examples are the so-called -Fp_i or type B 
models of the mode-coupling approach and disordered spin models with interac- 
tions among all possible p > 3-uples. These models also have only one additional 
time-scale. In these cases the x vs C plot, for C < qEA(7'), is a straight line 
of negative slope larger than -1/T [5], i.e. T < T^^^ < oo. The value of T^^) 
weakly depends on the working temperature and continuously increases from 
Tfi at the dynamic transition to a slightly higher value at T 0. Another in- 
teresting case in this group is a d-dimensional directed manifold embedded in 
an {N — )■ oo) -I- d-dimensional space under the effect of a random potential with 
short-range correlations. The susceptibility-correlation relations for different 
wave- vector dependent observables satisfy Xfc(C'fc) = Xfc=o(C'fc=o) = Xx=o{Cx=o) 
when times are such that they all evolve in the aging regime characterized by a 
single time scale (note that qff^ depends on k) and, thus, the effective tempera- 
ture is fc-independent [34]. 

Mean-field spin models with quenched disorder having a genuinely continuous 
second-order phase transition^ for example, the Sherrington-Kirkpatrick spin- 
glass or type A models of the mode-coupling approach. The random manifold 
problem with long-range potential correlations is of this type too. In these cases 
the dynamics has a continuity of timc-scalcs, ordered in an ultrametric fashion. 
For C < qEA{T) the x{C) plot is a non-trivial curve with local derivative larger 
than —1/T [6]. Each value of the effective temperature can be ascribed to 
a dynamic scale. The lowest value appears discontinuously (T*^^-* > Td) as one 
crosses Td and can be shown to decrease with decreasing temperature. The wave- 
vector independence holds for the Teff of the random manifold as well [34, 35] . 

Further support to the notion of effective temperatures comes from the study 
of the effect of quantum fluctuations on the same family of mean-field mod- 
els [36, 37]. The setting is one in which the system is in contact with a quantum 
environment at temperature T and the dynamics arc dissipative. Below a crit- 
ical surface (in the T, strength of quantum fluctuations, coupling to the bath 
phase diagram) that separates glassy from equilibrium phases, and in the slow 
dynamics regime, one finds a non-equilibrium relaxation with deviations from 
the quantum FDT, eq. (18). These are characterized by the replacement of the 
bath temperature by an effective temperature Teff. The effective temperature 
is again piecewise. It coincides with T when the symmetric correlation C is 
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larger than ^ea and Teff > T when C goes below ^ea- x{C) recovers, in the slow 
regime, the structure of the classical limit. This is a signal of a time-dependent 
decoherence effect. The slow modes' Teff depends on the working parameters 
and is higher than T even at T = 0. 

Numerical simulations, with Monte Carlo techniques or molecular dynamics 
of classical systems, have demonstrated that the classification can go beyond 
the mean- field limit. We shall discuss these tests in Sect. 5. 



3.4 Pre-asymptotic effects 

Expression (30) is asymptotic in time, a limit in which it is immaterial to con- 
struct the parametric plot by keeping fixed and letting t increase to infinity, 
or by keeping t fixed and letting t^j increase up to t. However, in numerical 
and experimental measurements one cannot reach this long-times limit and it is 
important to have the best possible control of pre-asymptotic effects. The most 
appropriate construction at finite times is the second option above [38] . Indeed, 
by using the fact that C{t,tw) is monotonic with respect to t^, one trades all 
ituS by Cs and keeps t as an independent variable. Variations with respect to 
tu, become variations with respect to C at i fixed, 

dt^ dc ^ ' n^^{t,c)dc rEFF(i,c) ' ^ ' 

and the slope of the pre-asymptotic x{t, C) plot at fixed t yields the searched 
pre-asymptotic Teff, as defined in cq. (20). Had we worked at fixed we 
would have failed to obtain Teff from the slope of the x(C, tw) plot. This is due 
to the asymmetry in the times involved in the definition (20). The difference 
disappears in the long tw and t limits with C varying in the interval [0, 1] or 
for stationary systems. In cases in which a normalization is needed, the choice 
of using C{t, t) as the normalization factor does not affect the slope of the plot 
constructed at t fixed. The FDRs discussed in Sect. 2.4 allow one to use this 
method with no extra computational effort in numerical simulations. 



3.5 Cooling rate dependence and finite size effects 

In the long-times dynamics of mean-ficld-like models Teff docs not depend per- 
manently upon the cooling procedure. The same models with finite number of 
degrees of freedom, and more refined ones beyond the mean-field approximation, 
should capture a cooling rate dependence that should also become manifest in 

Teff- 

The dynamics of models in the RFOT class approach, in the asymptotic limit 
t — >■ oo taken after A'^ — > cxd, a region of phase space, named the threshold [5] , that 
is higher than equilibrium in the free-energy landscape. Further decay is not 
possible in finite times with respect to N since diverging barriers separate the 
former from the latter. This is demonstrated by the fact that the asymptotic 
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values of averaged one-time intensive quantities, such as the energy density, 
take higher values than in equihbrium. Moreover, in between threshold and 
equilibrium a continuous set of metastable states also separated by diverging 
barriers exsit. Subsequent decay is achieved through activated processes in time- 
scales scaling with N, and the relaxation of one-time quantities is expected to 
cross over from power-law to logarithmic. 

Although the full dynamic solution in the activated regime has not been 
derived yet (it is too hard!), it is reasonable to imagine that the Teff values 
should be ordered with the higher on the threshold and the lower, T, in equilib- 
rium. The effective temperature should relax, in logarithmic time scales, from 
the threshold value to the bath temperature. Coming back to cooling rate de- 
pendencies, a slower cooling rate takes a finite size system below the threshold 
level and the deeper the slower the rate. Teff should follow this variation. These 
claims found some support in numerical simulations, see Sect. 5. 

3.5.1 An equilibrium interlude 

The replica theory of the statics of mean-field spin-glasses [30] necessitates the 
definition of a symmetric disorder-averaged functional order parameter, P{q), 
that measures the probability distribution of overlaps, Nq = X]i=i Si^ii between 
equilibrium configurations, {si\ and {ui}. The cumulative distribution and a 
further integral arc 



The out of equilibrium dynamics take place in a region of phase space that 
is different from the one where equilibrium states lie. Still, the classification 
of mean- field models according to rEFF(C') [or X{C)\ and x(C) coincides with 
the one arising from the replica analysis [30] of the number and organization of 
equilibrium states and their implications on x{q) and H(C) defined in eq. (32). 
Indeed, mean-field coarsening problems have two equilibrium states related by 
symmetry, 2P{q) = 6{\q\ - Qea) implying x{q) = if \q\ < Qea and x{q) = 1 if 
l^l > ^EA with (7ea the square of the order parameter (replica symmetric - RS - 
case) . Models with a RFOT are solved by a one step replica symmetry breaking 
(IRSB) Ansatz implying a proliferation of equilibrium states with special prop- 
erties. Two possibilities exist for any two equilibrium configurations {s;} and 
{<7i}: they may fall in the same or the reversed state and q ~ ±(7ea, respectively; 
or else they may fall in different states that turn out to be orthogonal and q = 0. 
The probability distribution is hence P{q) = xiS{q) -I- (1 — xi)/2d{\q\ — qE.\) and 
its integral yields x{q) = xi if \q\ < qsA and x{q) = 1 if \q\ > q^A- Finally, 
models of the SK type are solved by a full RSB Ansatz, have states with all 
kinds of overlaps, the P{q) has two delta contributions at i^EA and a symmet- 
ric continuous part in between. This functional form leads to an x{q) taking 
the value 1 for \q\ > (/ea and a continuous function of q for \q\ below Qea- In 





(32) 
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all mean-field cases the functional forms of equilibrium and out of equilibrium 
objects are similar, 

x{q) ^ X{C) , H(C7) o x{C) . (33) 

By this we mean that H has an FDT part in all cases and it is linear, with 
zero, constant and finite, and variable slope, in RS, IRSB-RFOT, and full RSB 
cases, respectively. In models of RFOT type the value of the breaking point ^ea 
and the parameter xi (the slope) are not the same dynamically and statically 
while in models of full RSB kind the coincidence is complete (up to a factor 2 
due to the global symmetry). The same applies to all higher moments of the 
equilibrium overlap q and out of equilibrium C distribution functions [5, 6, 35]. 
The coincidence for full RSB models was argued to apply beyond mean-field, in 
finite dimensional disordered spin systems, when the long-times limit is taken 
after the thermodynamic limit. Details of the reasoning, that is based on the 
assumptions of stochastic stability and the convergence of the out of equilibrium 
susceptibilities to the equilibrium can be found in [39] . 

These ideas were developed in spin models and one would like to extend 
them to atomic and molecular systems. However, overlaps in continuous particle 
models are difficult to define in a direct measurable way. Attempts based on 
weakly coupled real replicas were developed in [40]. This may allow one to 
extend the equilibrium •o- out of equilibrium connection to these systems as 
well. 

The relation between static H(C) and asymptotic out of equilibrium dynamic 
x(C) could apply in much more generality than previously suspected at the price 
of identifying finite time non-equilibrium, x(C, and finite size equilibrium, 
N(C, ^(^tu)), with the help of a time-dependent coherence length £,{t^) [41]. 

4 Requirements 

In this Section we list a number of conditions that the functional parameter Teff 
must satisfy to act as a temperature. 

4.1 Thermometer 

Any quantity to be defined as a non-equilibrium effective temperature must 
conform to the folklore. The first requirement is to be measurable with a ther- 
mometer weakly and statistically coupled to the system [3, 4]. This fact can 
be proven by studying the time-evolution of the thermometer coupled to M 
identical copies of the system, all of age t^, and evolving independently. The 
thermometer is ruled by a Langevin equation with a non-Markovian bath with 
statistics given by the system's correlation and response. It thus feels the system 
as a complex bath with its time scales, r'^')(i'), and temperatures, T^'^\ The T*^'^ 
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to be recorded is selected by tuning the internal time-scale of the thermometer 

to TW(i'). 

Such an experiment can be relatively easily realized in systems of particles in 
interaction, be them colloidal suspensions or powders. The thermometer can be 
a probe particle, the free and perturbed dynamics of which is followed in time. 
Diffusion is measured in free relaxation and mobility in the perturbed case. By 
comparing the two through an extended Einstein relation Teff of the medium, 
that is to say, the system of interest, is measured. Another possibility is to 
monitor the kinetic energy of the tracer (a quadratic variable) and associate it to 
the effective temperature of the environment (via equipartition). In both cases, 
by playing with the tracers' parameters, namely their mass, different regimes 
of relaxation are accessed and the T^^^s are measured. Consistency with the 
fact that the effective temperature should be an intensive variable requires the 
result to be independent of the shape of the tracers. These experiments have 
been performed numerically and experimentally yielding very good results in 
the former and somehow conflicting in the latter. Wc shall discuss them in 
Sect. 5. 

4.2 Role played by external baths 

The fluctuation dissipation ratio of an 'easy to equilibrate' system should ac- 
quire the temperature of its external environment. A rather simple though 
particularly illuminating problem that illustrates this idea is the non-Markov 
diffusion of a particle in a harmonic potential, simultaneously coupled to two 
baths, a fast one in equilibrium at temperature Tg"^ ~ T^"-' giving rise to white 
noise and instantaneous friction, and a slow one with an exponential memory 
kernel in equilibrium at temperature T^"^ = T^^\ This example can be taken as 
a schematic model for an internal degree of freedom in a slowly driven system 
and, in particular, the dynamics of a (possibly confined) Brownian particle in 
an out of equilibrium medium. If the slow bath is not stationary it can also 
be taken as a self-consistent equation for a variable in an aging system - exact 
for mean-field disordered spin models. The particle behaves as in equilibrium 
at r^°-' or T*-^-* depending on which bath it feels. The separation is made sharp 
by an adequate choice of the bath and potential energy parameters that pushes 
apart their own time scales. Similarly, an aging system with multiple effective 
temperatures becomes stationary in all time-scales with Tg"^ > T^ff and goes 
on aging in time-scales with T^*-* < Teff [4] . 

4.3 Observable dependencies 

In cases allowing for a thermodynamic interpretation, this intensive variable 
should be the same - partial equilibration - for all observables evolving in the 
same time-scale and interacting strongly enough. A concrete check of this fea- 
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ture was performed within solvable glassy models. Two glasses in contact with 
a bath in equilibrium at temperature T, each of them with a piecewise rEFF(C') 
of the form 



rrifSYST 1,2) _ ) T if C > ^ea' 



1.2) 
;a 
,(1-2) 



with r(i) ^ r(2), were chosen. The experiment of setting two observables in 
contact is reproduced by introducing a small linear coupling between micro- 
scopic variables of the two systems. Above a critical (though small) value of 
the coupling strength the systems arrange their time-scales so as to partially 
come to equilibrium and the effective temperatures below ^ea equalize. Below 
the critical strength T^^'^^ remain unaltered [4]. 

The observable independence of the fluctuation dissipation ratio was inves- 
tigated in a variety of out of equilibrium situations but no clear rationale as 
to when this holds was found yet. An intriguing proposal was put forward by 
Martens et al. who argued that the observable independence and hence the 
interpretation in terms of Teff is related to the uniformity of the phase space 
pdf on the hyper-surface of constant energy reached dynamically [21]. These 
authors validated this idea in a few simple toy models relaxing in a single time- 
scale (glassy systems excluded) . Consistently, this proposal applies in mean field 
disordered models. Moreover, they proposed that the observable dependence is 
proportional to the square root of the difference between the Shannon entropy 
of the dynamic state and the equilibrium one, a conjecture that deserves further 
investigation. 



4.4 Intuitive properties 

One may wonder why in all studied un-driven systems there is a two-time regime 
in which C decays from its equal times value to ^ea and FDT holds. This 
problem admits a physical explanation - thermal fluctuations within domains, 
rapid vibrations within cages, etc. - that we shall expose in Sect. 5, a free- 
energy landscape explanation that we shall discuss in Sect. 4.8 but also a formal 
explanation: there exists a bound - on the difference between left and right hand 
sides of eq. (13) - that vanishes in the first regime of relaxation [42]. 

An intuitive property of Teff is that it loosely represents the disorder level of 
the system. This idea translates into Teff being higher or lower than the working 
temperature T when the initial state is equilibrium in the disordered phase or 
at a lower temperature than the quenching value. These features are realized in 
all cases studied so far, mean- field or finite dimensional alike [43, 44, 45]. 

This slated property matches the 'fictive temperature', Tf, ideas that date 
back to the 40s at least [8] and have developed ever since, as explained in [9]. 
Tf is usually introduced by assuming that when a liquid falls out of equilibrium 
in the glass transition region its structure gets 'frozen' at a Tf that is higher 
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than r, depends upon the coohng rate and in particular on T, and deep below 
the transition range approaches Tg. Several definitions in terms of the enthalpy 
or the thermal expansion coefficient have been given and they do not necessar- 
ily coincide. The Active temperature is hence a phenomenological convenience 
and acts essentially as a parameter in an out-of-cquilibrium 'equation of state'. 
In contrast, the effective temperature is defined in terms of fluctuations and 
responses, it can be measured directly, and plays a role that is closer to the 
thermodynamical one (although in not all possible out of equilibrium systems 
but in a class yet to be defined precisely). 

4.5 Fluctuation theorem 

The fluctuation theorem concerns the fluctuations of the entropy production 
rate in the stationary non-equilibrium state of a driven dynamical system [46] . 
It applies to systems that approach equilibrium when the forcing is switched off. 
Of interest in the context of glassy systems is to know the fate of the fluctuation 
theorem if the system evolves out of equilibrium even in the absence of the 
external drive, and whether Teff enters its modified version. These questions 
were addressed by Sellitto who analyzed the fiuctuations of entropy production 
in a kinctically constrained lattice model [47] and by Crisanti and Ritort who 
discussed the interplay between Teff and a fluctuation theorem on heat-exchange 
between the system and the environment in the random-orthogonal model (a 
case in the RFOT class) [48]. From a different perspective, closer to the one 
in [15], the fluctuation theorem and Teff were studied in [49] within two related 
problems: mean-flcld glassy models and a Langevin process with a number 
of equilibrium thermal baths with different time-scales and temperatures as 
the one discussed in Sec. 4.2. Firstly, it was shown that the work done at 
frequency uj by conservative and non-conservative forces is weighted by the 
effective temperature (instead of the temperature of the bath) at the same 
frequency. The work of the conservative forces produces entropy if the bath 
is out of equilibrium since the nonlinear interaction couples modes at different 
frequency which arc at different temperature, thus producing an energy flow 
between these modes. Secondly, it was proven that the entropy production rate 
satisfies a fluctuation theorem. Thirdly, extensions of the Green-Kubo relations 
for transport coefficients were derived. Fourthly, a feasible way to measure 
Teff by exploiting the modifled fluctuation theorem was discussed. There is no 
satisfactory numerical test of these ideas in non-mean-field glassy systems yet 
but a preliminary study will be discussed in Sect. 5.7. 

4.6 Non-linear effects 

In a series of papers Hayashi, Sasa et al. [50] investigated the notion of an 
effective temperature in classical non-equilibrium steady state (NESS) focusing 
on a strongly perturbed Id white-noise Langevin system in which a Brownian 
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particle is subject to a spatially constant driving force / and a periodic potential 
U{x). These authors extended the definition of Teff that uses the Einstein 
relation between diffusion coefficient and differential mobility beyond the linear 
response regime, Teff(/) = D{f)/fj,{f). They addressed the question of the 
physical significance of Teff in at least three different ways: by showing that it 
plays the role of a temperature in a large-scale description; by proving that such 
an out-of-equilibrium system, used as a thermostat for a Hamiltonian system, 
is able to transfer its effective temperature as kinetic energy; and by conducting 
a heat conduction experiment. The results of these tests are consistent with a 
thermodynamic interpretation. 

4.7 Local measurements 

The fluctuations-dissipation theorem relates the averaged local response func- 
tion and fluctuations measured in equilibrium at any spatial scale. What should 
local measurements yield out of equilibrium? 

4.7.1 Quenched disorder induced fluctuations 

In systems with quenched random interactions spatial fluctuations in noise- 
averaged quantities arc dictated by the local disorder. This fact has been 
known for decades; for instance, the fast and slow character of spins in dis- 
ordered magnets induced by their quenched environment, give rise to Griffiths 
phenomena (free-energy singularities, phases and slow relaxation). In the con- 
text of glassy dynamics and Teff, Montanari and Ricci-Tersenghi [51] showed 
that in disordered spin models defined on random graphs spins of two types 
exist: paramagnetic and glassy ones with the former following fast equilibrium 
dynamics and the latter having a non-trivial relaxation peculiar to their envi- 
ronment. Strictly local deviations from FDT were characterized by an effective 
temperature that can be obtained with a replica calculation along the lines dis- 
cussed in Sect. 3.5.1. Pretty convincing arguments - although not completely 
rigorous - establish that in these models the strictly local effective temperature 
must diffuse to become site independent 

= (34) 

TEFF(a) dc, 

Numerical studies of this kind of (absence of) Teff fluctuations in the 3d EA 
model [52] will be presented in Sect. 5. 

4.7.2 Noise induced fluctuations 

In a series of papers Chamon et al. proposed a theoretical framework based 
on global time-reparametrization invariance that explains the origin of dynamic 
fluctuations in generic - not necessarily quenched disordered - glassy systems 
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with a separation of time-scales of the kind explained in Sect. 3.2 [53, 54, 55]. 
Such type of invariance had been known to exist in the out of equilibrium dy- 
namics of mean-field disordered models and it was later shown to carry through 
to the causal asymptotic dynamics of finite d infinite size spin-glasses, under 
the assumption of a slow dynamics with a separation of time-scales. The ap- 
proach reviewed in [56] focuses on the fiuctuations induced by the noise at 
coarse-grained length scales. The invariance acquires a physical meaning and 
it implies that one can easily change the clock h{t) characterizing the scaling 
of the global correlation and linear response by applying infinitely weak per- 
turbations that couple to the zero mode, or with a noise-induced fiuctuation. 
An illustration of this property is the fact that the aging relaxation dynamics 
of glassy systems is rendered stationary by a weak perturbing force that does 
not derive from a potential while the xi^) relation in the slow regime is not 
much modified [57]. The same argument, applied to the fluctuations, implies 
that easy fluctuations should be realized as local changes in time, t — > hr{t), 

MI ^ 1 + - 1') ^ 1 , i') 

hr{t') ~ dt'lnhrit') Tr{t') ^ ' 

that intervene in the local correlation and integrated linear response. Age mea- 
sures fluctuate from point to point with younger and older pieces coexisting 
at the same values of the two laboratory times. The two-time coarse-grained 
obscrvables, i.e. Cr{t,tm) and Xrit^tw), have a slow and a fast contribution, 
the former characterized by scaling functions fc and f^^ that act as 'massive 
variables' in the sense that they are not expected to fluctuate in the scaling 
limit 5 <C £ ^ with 6 the microscopic length-scale (lattice-spacing or 

inter- particle distance), i the coarse- graining length (to be chosen), and ^ the 
dynamically generated correlation length. Within this picture the parametric 
construction Xr(C'r) falls on the master curve for the global quantities but could 
be advanced or retarded with respect to the global value with a uniform effective 
temperature at fixed Cr value: 

T'EFFr(C'r) = T^FviCr) • ("^6) 

This equality is non-trivial when the effective temperature is finite and it con- 
forms to the concept that the effective temperature for different obscrvables 
(the different regions in the sample in this case) must equalize if they evolve in 
the same time-scale {i.e. take the same value of Cr). Timc-rcparametrization 
invariance is not expected to hold in cases in which Teff diverges [58] , notably, 
phase ordering kinetics (see Sect. 5.2). 

Equation (36) with T^FFiCr) finite implies that the two-time variances of 
composite fields, the averages of which yield the local correlation and linear 
response should have the same scaling with times [59]. These objects are easier 
to measure than the joint probability distribution function (pdf) of Xr and Cr- 
Numerical results shall be discussed in Sect. 5. 
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4.8 Landscapes and thermodynamics 

Complex systems' dynamics are sometimes interpreted as the wandering of a 
representative point in a phase space endowed with a comphcated free-energy 
density landscape. The existence of an equilibrium-like relaxation at short- 
time differences suggests the distinction between 'transverse' and 'longitudinal' 
directions in the landscape, with the former being confining and close to, say, 
'harmonic' and the latter giving rise to the slow out of equilibrium structural 
relaxation. Questions are posed as to what are these directions, which is the 
dynamic process taking the system along these directions - diffusive, activated, 
etc. — which is their volume in phase space, and so on. 

The appearance of an effective temperature suggests some form of ergodocity 
and it becomes tempting to relate Teff to some kind of microcanonic tempera- 
ture defined from the volume of phase space visited during the out of equilibrium 
excursion. The construction of a thermodynamics in which Teff played a role 
is the natural step to follow. In this Section we describe the successful devel- 
opment of this program in models of the RFOT sort and how this construction 
remains a suggestive picture in finite dimensional cases. 

4.8.1 Thermal systems 

Describing super-cooled liquids and glasses in terms of potential energy land- 
scapes dates back to Goldstein who proposed to think of a system's trajec- 
tory in phase space as a succession of steps among potential energy basins [60] . 
This idea developed into the inherent structure (IS) statistical mechanics frame- 
work of Stillinger and collaborators [61]. In this approach each configuration is 
mapped onto a local minimum of the potential energy through a minimization 
process implemented, for example, by a quench to T = (steepest descent). 
The inherent structure is, then, the configuration reached asymptotically and 
all configurations flowing to it constitute its basin of attraction. Although a 
priori simple, this proposal hides a number of ambiguities such as the fact that 
the ISs depend on the microscopic dynamics {e.g. single spin flip vs. cluster 
spin flip in a spin system), some decision making is needed in cases in which 
the T = dynamics could follow different directions, etc. The proposal is to 
re-order and approximate the partition function as a sum over IS energy levels 
(including their degeneracy) times a /3-dependent factor with all contributions 
from the rest of the configurations - associated to vibrations or the fast relax- 
ation - basically describing the free-energy of the liquid/glass constrained to 
one typical basin. The assumption is that the IS partition function describes 
the thermodynamic properties of the state reached at very long times. 

Independently, Thouless, Anderson and Palmer (TAP) [31] and de Dominicis 
and Young [62] showed that the equilibrium properties of fully-connected spin 
disordered models can be described with a local order-parameter dependent 
free-energy landscape. Averaged observables in equilibrium are expressed as 
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weighted sums over the free-energies of the TAP free-energy saddle-points: 



in agreement with resuUs from rephca and cavity methods. The sums in the 
second member run over aU stationary points of the TAP free-energy landscape. 
They transform into integrals over free-energies at the price of introducing the 
number of stationary points at given /, with iVE(/) = hiAf{f) the complexity, 
or configurational entropy at free-energy density /. This construction is exact 
for mean-field models. In the N ^ oo limit the integral is evaluated by saddle- 
point: 

1 



T df 



(38) 



In RFOT models the disordered state dominates above and /eq ~ /p^; in 
between and states that are not minima of the TAP free-energy control 
the integral since their number is sufficiently large, /eq = fsp — TS(/sp); finally, 
at Ts the configurational entropy vanishes (entropy crisis) and the lowest lying, 
now glassy, states dominate. 

The microcanonical vision of the effective temperature suggests to check 
whether below 

1 _ as(/) 



Teff Of 



(39) 

/th 



where Teff is the value of the slow modes effective temperature and /th is the 
free-energy at the threshold, the level reached dynamically by an infinite system 
after a quench from high T. This is indeed the case in RFOT type models in 
the thermodynamic limit. Systems with large but finite size relax below the 
threshold and slowly approach equilibrium. Nieuwenhuizen conjectured that in 
these cases both Teff and E(/) acquire a time-dependence in such a way that 
eq. (39) remains valid with S(/, <) the complexity of the TAP states that are 
relevant at time t [63]. Moreover, in finite-dimensional systems with short-range 
interactions, barrier heights and lifetimes are finite at finite temperature and 
metastability becomes a matter of time scales. A recipe to compute S(/, t) in 
these cases was given in [64]. 

Once a separation of modes into fast (vibrational) and slow (structural) is 
made in either an approximate (in finite dimensional models with short-range 
interactions) or exact (in mean-field cases) way, thermodynamic potentials that 
involve Teff can be easily constructed and thermodynamic relations derived. 
This has been done in the IS formalism [61] and in a framework that is closer 
to the TAP one [63]. 

The inadequacy of the IS approach to describe the dynamics of coarsen- 
ing and kinetically constrained models has been explained in [65]. Its limits of 
applicability in molecular glasses were also discussed. Nevertheless, since it is 
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not evident how to access a free-energy landscape concretely, numerical efforts 
have focused on the characterization of the potential energy landscape. While 
increased computational facilities gave access to an exhaustive enumeration of 
ISs in small clusters and proteins, the calculations remain incomplete for macro- 
scopic systems. A connection between Teff and the IS complexity, that has to 
be taken with the caveats mentioned above, was discussed in [66], see Sect. 5. 



4.8.2 Athermal systems: Edwards ensemble 

The constituents of thermal systems exchange energy with the components of 
their environment and this exchange has an effect on their motion. The con- 
stituents of athermal systems are much larger than the ones of their surroundings 
and the energy received from the bath is irrelevant. Dissipation occurs via en- 
ergy flow from the particles to internal degrees of freedom that are excluded 
from the description. Granular matter is the prototype. 

In spite of the very different microscopic dynamics, the meso/macroscopic 
dynamics of gently perturbed dense granular matter share many points in com- 
mon with the ones of more conventional glassy systems. A detailed description 
of the experiments demonstrating these facts is given in [67]. 

In search for a statistical mechanics description of these systems, Edwards [2] 
proposed to use the volume, V, as the macroscopic conserved quantity, and the 
blocked states - defined as those in which every particle is unable to move - as 
the set of relevant equiprobable configurations, P{C, V) — {V)d{C)6{V{C) — 
V)), where 0{C) is an indicator function that equals one if the configuration is 
blocked and zero otherwise, V(C) is the volume as a function of the configuration, 
and ^1{V) = J dC P{C, V) is the volume in configuration space occupied by the 
blocked states. An entropy and compactivity are next defined as 

S{V) = -J2PiC,V)\nPiC,V)^\nn{V), J— ^ ^ (40) 

g ^ EDW O r 

respectively. A strong hypothesis in this description is that all blocked con- 
figurations are treated on equal footing and that any distinction of dynamic 
origin is disregarded. Moreover, the number of conserved macroscopic variables 
needed to correctly describe the system is not obvious a priori. If the grain 
configurations arc also characterized by an energy one should then enlarge the 
description and define 

P{C,V, E) = n-\V, E)e{C)6{V{C) - V))S{H{C) - E)) , 

SiV,E) = -J2PiC,V,E)\nPiC,V,E), 1 ^ dS{VE) ^^^^ 

^ -L EDW O-Cj 

The similarity between the entropy of blocked states and the zero-tem- 
perature complexity of glass theory is rather obvious. Indeed, the first successful 
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check of Edwards' hypothesis was achieved in RFOT models at T ~ [7]. In 
these cases the energy is the relevant macroscopic variable. One identifies all 
energy minima (the blocked configurations in a gradient descent dynamics), cal- 
culates I/Tedw and shows that it coincides with Teff = r'^-* in the slow aging 
regime (C < ^ea) [5]. Moreover, all intensive observables are given by their flat 
average over the threshold level Cth = linit_).oo li™-Ar-!.oo e(i). At finite T the 
connection can be extended at the expense of using the free-energy instead of 
the energy, as explained in Sect. 4.8.1. 

Numerical tests in finite dimensional kinetically constrained lattice gases [68] , 
microscopic models of sheared granular matter including some of the subtleties 
of frictional forces [69], spin-glasses with athermal driving between blocked 
states [70] and a particle deposition model [71], gave positive results. In all 
these cases Edwards measure is able to correctly reproduce the sampling of the 
phase space generated by the out of equilibrium dynamics. Nevertheless, this 
description docs not apply to every problem with some kind of slow dynamics. 
In [68] the counter-example is the domain growth of a 3d random field Ising 
model, a case in which the properties of a long-time configuration of (low) en- 
ergy is not well reproduced by the typical blocked, by domain-wall pinning by 
disorder, configuration of the same energy. In [72] , instead, the analytically solv- 
able one-spin flip dynamics of the Id Ising chain is used to display quantitative 
and qualitative discrepancies between the dynamic treatment and the averaging 
over an a priori probability measure of Edwards type (and refinements). At 
the mean-field level the SK model and the like do not admit a simple relation 
between configurational entropy and Teff either. 

A careful account of the experimental subtleties involved in trying to put 
Edwards' hypothesis to the test, and eventually verifying whether Tedw — Teff, 
is given in [67]. The question remains open especially due to the difficulty in 
identifying the relevant extensive and intensive thermodynamic parameters [73] . 
All in all, the approach, very close to the IS and TAP constructions, is intriguing 
although not justified from first principles yet and its limits of validity remain 
to be set. 

5 Measurements 

In this Section we discuss measurements of FDT violations and tests of the 
effective temperature notion in a variety of physical systems out of equilibrium. 
Since we cannot make the description exhaustive we simply select a number 
of representative cases that we hope will give a correct idea of the level of 
development reached in the field. 
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5.1 Diffusion 



The dynamics of a particle in a potential and subject to a complex environ- 
ment (colored noise or baths with several time-scales and temperatures) has a 
pedagogical interest but also admits an experimental realization in the form of 
Brownian particles immersed in, e.g., colloidal suspensions and controlled by 
optical tweezers. 

A particle coupled to a bath in equilibrium at temperature T with noise- 
noise correlations of type (^(t)C(t')) c>c (i — t')~"-~^ , < a < 2, and under no 
external forces, performs normal or anomalous diffusion depending on a. The 
fluctuation-dissipation ratio, cq. (20), for t>t'\s [74] 

, TR^^{t,t') _ D{t~t') 
^"^^'^ ^ " d,CUt,t') " D{t-i') + D{t') ' ^^^> 

with the diffusion coefficient D{t) = 1/2 d{x^{t))/dt ~ for a ^ 1 and D(t) = 
ct for a = 1. In the colored noise cases X^^ is a non-trivial function of times 
and it does not seem to admit a thermodynamic interpretation. Still, for later 
reference we consider the long times limit: 

subOhmic, 
Ohmic, 
superOhmic. 



lim lim = X^^{t,t') ~ < 1/2 a = 1 Ohmic 






a 


< 1 


1/2 


a 


= 1 


1 


a 


> 1 



Another illustrative example is the non-Markovian diffusion of a particle 
in a harmonic potential and subject to different external baths [4, 49, 75]. As 
already explained in Sect. 4.2 this simple system allows one to show how different 
environments can impose their temperatures on different dynamic regimes felt 
by the particle. Tests of other definitions of out of equilibrium temperatures 
in this simple case confirmed that the definition that appears to have the most 
sensible behaviour is the one stemming from the long-time limit of the relations 
between induced and spontaneous fluctuations [75]. All other definitions yield 
results that are more difficult to rationalize: in most cases one simply finds the 
temperature of the fast bath and in some cases, as with the static limit in [76], 
one incorrectly mixes different time regimes even when their time-scales are well 
separated. 



5.2 Coarsening 

When a system is taken across a second order phase transition into an ordered 
phase with, say, two equilibrium states related by symmetry, it tends to or- 
der locally in each of the two but, globally, it remains disordered. As time 
elapses the ordered regions grow and the system reaches a scaling regime in 
which time-dependencies enter only through a typical growing length, L{t). Fi- 
nite dimensional coarsening systems have been studied in great detail from the 
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effective temperature perspective (see [13]). In this context, it is imperative to 
distinguish cases with a finite temperature phase transition and spontaneous 
symmetry breaking from those with ordered equihbrium at T = only. Some 
representative examples of the former are the clean or dirty 2d Ising model with 
conserved and non-conserved order parameter. An instance of the latter is the 
Glauber Ising chain and we postpone its discussion to Sect. 5.4. 

Let us focus on scalar systems with discrete broken symmetry. When time- 
differences are short with respect to the typical growing length L{t^), domain 
walls remain basically static and the only variation is due to thermal fluctua- 
tions on the walls and, more importantly, within the domains. This regime is 
stationary, and induced and spontaneous fluctuations are linked by the FDT. 
At longer time-differences domain walls move and observables display the out 
of equilibrium character of the system. 

The correlation and total susceptibility in the iu, — )■ oo limit separate in two 
contributions C{t,U,) = C'^it - t„) + C^'^\t,t^) and xit,tw) = x''^(i - t^') + 
x'^"'^' (t, ^u,). Numerical studies of Teff focused on the parametric construction 
x{C,tyj) at fixed and finite where the chosen observable is the spin itself. 
The resulting plot has a linear piece with slope — as in eq. (13), that goes 
below C = Qea = rn^ and, consistently, beyond X = [1 ~ 'm?]/T. The addi- 
tional equilibrium contribution is due to the equilibrium response of the domain 
walls that exist with finite density at any finite tw In the truly asymptotic 
limit their density vanishes and their contribution disappears. Consequently, 
limj^^oo x(C', ^iii) = C^"^ > Qea satisfies FDT and it is entirely due to fluctua- 
tions within the domains. In cases with L{t) ~ t^/^"^, the slow terms take the 
scaling forms 



defects p{t) ~ L(t)~^ ~ i-n/z^ with ri = 1 for scalar and n = 2 for vector order 
parameter. Although this seems plausible is instead d-dependent. Another 
conjecture is [13] 



(11 is the dimension at which vanishes and may coincide with the lower 
critical dimension. One finds ^ I'm the Ising model, = 1 in the Gaussian 
approximation of Ohta, Jasnow and Kawasaki, and = 2 in the 0{N) model 
in the large N limit, du is the dimension at which a^^ becomes d-indcpendcnt 
and it does not necessarily coincide with the upper critical dimension. One finds 
dc/ = 3 in the Ising model, dc/ = 2 in the Gaussian approximation, and du = A 
in the large N 0{N) model. It was then suggested that djj might be the highest 
d at which interfaces roughen. In all cases in which > 0, T^ff ^ oo- This 





d < du , 
d = du , 
d > du ■ 



(44) 
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result was confirmed by studies of second order FDRs in tlie 2d Ising model 
that showed the existence of stationary contributions verifying the non-linear 
equilibrium relation and aging terms that satisfy scaling and yield Teff — > oo as 
in the linear case [25]. The approach by Henkel et al. based on the conjecture 
that the response function transforms covariantly under the group of local scale 
transformations, fixes the form of the scaling function but not the exponent 
a-)^ [77] and does not make predictions on Teff- The coincidence between statics 
and dynamics, see Sect. 3.5.1, holds in these cases [13]. 

Noise induced spatial fluctuations in the effective temperature of clean coars- 
ening systems were analyzed in the large N 0{N) model [58] and with numerical 
simulations [78] . The first study shows that time-reparametrization invariance 
is not realized and that Teff is trivially non-fiuctuating in this quasi-quadratic 
model. The second analysis presents a conjecture on the behaviour of the av- 
erage over local (coarse-grained) susceptibility at fixed local (coarse-grained) 
correlation that consistently vanishes in coarsening (but is more interesting in 
critical dynamics as we shall discuss in Sect. 5.3). 

The results gathered so far and summarized in the conjecture (44) imply 
that the FD ratio vanishes and thus Teff diverges in quenches into the ordered 
phase of systems above their lower critical dimension. This 

5.3 Critical dynamics 

The non-equilibrium dynamics following a quench from the disordered state to 
the critical point consists in the growth of the dynamical correlation length, 
^(t) ~ This length docs not characterize the size of well defined domains 

but the size of a self-similar structure of domains within domains, typical of 
equilibrium at the critical point. A continuum of finite time-scales associated 
to different wave- vectors, t^''^ ~ k~''"', exists with only the fc — ?► diverging. 
At any finite time t, critical fluctuations of large wave- vectors, k^{t) ^ 1, are in 
almost equilibrium, while those with small wave- vectors, k^(t) <IC 1, retain the 
non-equilibrium character of the initial condition. This finite-time separation, 
and the fact that the order parameter vanishes, leads to the multiplicative scaling 
forms 

c(t,i^) ~e(^-^^„ )-'+'-" .fc[m/at^),^o/atra)] , 

x{t,t^)^p-at~t,,)-'^'-^ /xK(^)/e(^»),^o/e(^^.)] , 

with the microscopic length ensuring the normalization of the correlation 
and the fact that x vanishes at equal times. These forms imply that beyond 
the initial equilibrium part, the x(C') plot assumes a non-trivial shape that, 
however, progressively disappears and approaches the equilibrium linear form at 
all C > 0. The limit C = is distinct and the limiting parameter X°° should be 
non-trivial and universal in the sense of the renormalization group [79] . Whether 
this one can be interpreted as a temperature is a different issue that has been 
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only partially discussed. For this reason, we keep the notation X°° (instead of 
Teff) in most of this section. 

The correct estimation of X°° has to take into account that the number 
of out of equilibrium modes decreases in the course of time (contrary to what 
happens in the random manifold problem in the large A'' limit, for example) . The 
best determination of X°° is achieved by selecting the fc — >■ mode. A thorough 
review of the universality properties of X°° found with the perturbative field- 
theoretical approach and some exact solutions to simple models, as well as the 
comparison to numerical estimates, is given in [12]. In conclusion, X°° is a 
universal quantity that does not depend on the observable - as checked for a 
large family of them - but recalls certain features of the initial condition [45] 
and the correlations of the environment [80]. In the scalar model one finds the 
diffusive results, eq. (43), at the Gaussian level and corrections when higher 
orders are taken into account. For example X°° = 0.30(5) in d = 2, X°° = 
0.429(6) in d = 3 for a quench from a disordered state, white noise and up to 
second order inA-d. Instead, ~ 0.78 {d = 3) and = 0.75 {d = 2) if the 
initial state is magnetized. A larger X°° implies a lower = Tc/X°° and the 
comparison between these values conforms to the intuitive idea that an ordered 
initial state leads to a lower effective temperature than a disordered one. 

A different type of critical phenomena (infinite order) arises in the 2c? XY 
model. The magnetic order parameter vanishes at all T but there is a low-T 
critical phase with quasi long-range order (power-law decaying spatial corre- 
lations) that is destroyed at Tkt where vortices proliferate and restore a fi- 
nite correlation length. Out of equilibrium the critical scaling forms apply al- 
though with a temperature-dependent exponent, ri{T), and a growing length 
scale ^(t) ~ (t/lnt)^/^ (the logarithm is the effect of vortices). The role of 
the EA order parameter is played by the asymptotically vanishing function 
(t^/lnt„)-''(^)/2 and the crossover between equilibrium and out of equilibrium 
regimes takes place at a tu,-dcpendcnt value of the correlation. The x(C', t^) plot 
at finite is curved, it does not reach a non-trivial master curve for — > oo, 
but TEFF{t,t.ui) = fx[^{t)/^{tw)]. Quenches from the disordered phase, Tq > Tkt 
and heating from a To = ground state to T < Tkt demonstrate that the slow 
modes' Teff depends on the initial state and it is higher (lower) when Tg > T 
(To < T) [43]. We allow ourselves to use Teff in this case since these results 
point in the direction of justifying its thermodynamic meaning. Similar results 
were obtained for 1 + 1 elastic manifolds with and without quenched disorder 
(see [44] and refs. therein). As the dynamic-static link is concerned, Berthier 
et al. evinced that the extension to finite-times finite-sizes works, at least at 
not too high Ts where free vortices inherited from the initial condition arc still 
present. The coexistence of a single time scale in the aging regime together with 
a smooth and time-dependent x(C, i^,) plot arises naturally in a critical regime 
and it is due to the lack of sharp time-scale separation. 

The exact calculation of the joint probability distribution of the finite- 
size correlation and linear response in the spherical ferromagnet quenched to 
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its critical temperature was given in [81]. The results prove that these fluc- 
tuations are not linked in a manner akin to the relation between the aver- 
aged quantities, as proposed in [56], see Sect. 4.7.2, for glassy dynamics. The 
correlation-susceptibility fluctuations in non-disordered finite-dimensional fer- 
romagnets quenched to the critical point were examined in [78] where it was 
shown that the restricted average of the susceptibility, at fixed value of the two- 
time overlap between system configurations, obeys a scaling form. Within the 
numerical accuracy the slope of the scaling function yields, in the asymptotic 
limit of mostly separated times, the universal value X°° . 

The first experiments testing fluctuation dissipation deviations in a liquid 
crystal quenched to its critical point appeared recently and the results are fully 
consistent with what has been discussed above [82]. 

Although many evaluations of in a myriad of models tend to confirm 
that it mostly behaves as a critical property [12], the thermodynamic nature 
of this parameter has not been explored in full extent yet. Measurements with 
thermometers and connections to microcanonical definitions have not been per- 
formed at critical points. 

5.4 Quenches to the lower critical dimension 

The kinetic Glauber-Ising spin chain is the protoype of a dynamic model at its 
lower critical dimension. Taking advantage of the fact that this is one of the very 
few exactly solvable models of non-equilibrium statistical mechanics, several 
issues concerning the effective temperature interpretation have been addressed 
in this case, as the observable dependencies. 

After a quench from Tq — > oo to T = the factor Xs{t, t'), associated to the 
spin correlation and susceptibility, is < 1 and its value in the limit Cg ^ 
evolves smoothly from 1/2 (as in models characterized by simple diffusion such 
as the random walk or the Gaussian model [17]) to 1 (equilibrium) as t/req grows 
from to oo [l/re^ = 1 — tanh(2 J/T) is the smallest eigenvalue of the master 
equation operator]. Moreover, X^ is an exclusive function of the auto-correlation 
Cs as in more complex instances of glassy behaviour [84] . 

The value for the long-wavelength analogue, the fluctuating magnetization, 
X^, is identical to the local value Xg- The physical origin of the local-global 
correspondence, which can also be obtained by field-theoretic arguments [12], is 
that the long wavelength Fourier components dominate the long-time behaviour 
of both quantities. In contrast, obscrvablcs that arc sensitive to the domain 
wall motion have X^ = [83] , the difference residing on the interplay between 
criticality and coarsening, a peculiar feature of models with Tc = [13, 83]. 

The dependence on the initial condition is also interesting. A non-zero initial 
magnetization does not change the value of X^ at T = 0. Instead, demagne- 
tized initial conditions with strong correlations between spins so that only a 
finite number of domain walls exist in the system, yield X"^ = (the same re- 
sult is found in the spherical fcrromagnct) [38]. The deviations from non-linear 
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FDTs have not been fully analyzed yet. 

The static-dynamics connection [39] sketched in Sect. 3.5.1 does not hold 
in the Id Ising chain [13] and the non-trivial x(C') cannot be used to infer the 
properties of the equilibrium state. Indeed, the aging part of the response is 
finite asymptotically while the equilibrium P{q) has a double-delta (RS) struc- 
ture as in higher dimensions. The reason for the failure is that the hypotheses 
used to derive the connection are not fulfilled. 

The large N 0{N) model m D = 2 shares many common features with 
the phenomenology described above [58] although it has not been studied in as 
much detail. 

To sum up, a quench to T = at the lower critical dimension does not seem 
to be the dimensional continuation of a line of critical quenches in the (T, d) 
plane (as often implicitly assumed), but the continuation of a line of T = 
quenches: the system behaves as in the coarsening regime, although X°° ^ 
for observablcs that do not focus on the domain wall dynamics [13]. 

5.5 Relaxation in structural glasses 

In particle glassy systems a separation of time-scales exists although it is not 
as sharp as in mean-field models or coarsening systems, at least within sim- 
ulational and experimental time-scales. In atomic glasses the existence of an 
FDT part implies that the rapid particle vibrations within the cages occur in 
equilibrium while the structural relaxation is of a different out of equilibrium 
kind, and it is not necessarily ruled by the temperature of the bath. Tests of the 
thermodynamic origin of fluctuation-dissipation violations in the aging regime 
of these systems were carried through in much more detail and we summarize 
them below. 

5.5.1 Simulations of microscopic models 

Mono-atomic and binary Lcnnard- Jones mixtures, soft sphere systems, and the 
BKS potential for silica are standard models for glass forming liquids. Both 
Monte Carlo and molecular dynamics simulations [85, 86] suggest that the 
three first cases belong to the RFOT class of systems defined in Sect. 3.3 with 
constant in the aging regime. T*^^) depends weakly on the bath 
temperature and systems' parameters but it does not on the preparation proto- 
col as demonstrated by measurements after quenches and crunches [86] or the 
microscopic dynamics [87]. Tests of partial equilibration between fluctuations 
at different wave- vector gave positive results [85]. Importantly enough, these 
models have a well deflned equilibrium behaviour and their energy density is 
naturally bounded. Of special interest is the numerical method devised to com- 
pute linear responses in molecular systems with high precision that allowed one 
to resolve the paradoxical behavior previously reported for silica [88] . 
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Numerical evidence for a slow decrease in time of the configurational temper- 
ature, as defined in eq. (39), although with the inherent structure complexity, 
is in agreement with the idea of the system's representative point penetrating 
below the threshold in the (frce)-energy landscape [66]. 

The ratchet effect of an asymmetric intruder in an aging glass was studied 
numerically in [89]. The energy flowing from slow to fast modes is rectified to 
produce directed motion. The (sub) velocity of the intruder grows monotonically 
with Teff/T" and this current could be used to measure Teff- 

5.5.2 Kinetically constrained models 

Kinetically constrained models are toy models of the glassy phenomenon [90, 
91]. Their equilibrium measure is just the Boltzmann factor of independent 
variables and correlations only reflect the hard core constraint. Still, many 
dynamic properties of glass forming liquids and glasses arc captured by these 
models, due to the sluggishness introduced by the constrained dynamic rules. 
The literature on kinetically constrained models is vast; a recent review with 
tests of Teff is [91]. In short, non-monotonic low-temperature response functions 
were initially taken as evidence against the existence of effective temperatures 
in these systems. The confusion aroscd from the incorrect construction of the 
x(C) plot by using t^^ instead of t fixed (see Sect. 3.4) that led to the incorrect 
treatment of the transient regime. Still, even this taken into account, a large 
number of observables have negative fluctuation-dissipation ratios; this might 
be related to the fact that these models do not have a proper thermodynamics. 

5.5.3 Experiments 

Grigera and Isracloff were the first to measure FDT violations in glasses by 
comparing dielectric susceptibility and polarization noise in glycerol at T = 
179. 8K, i.e. relatively close to Tg ~ 196K [92]. At fixed measuring frequency 
uj ~ 8Hz, they found an effective temperature that slowly diminishes from 
Teff — 185K to roughly 180K in 10^ sec, that is to say in the order of days! 
This pioneering experiment in such a traditional glass former has not had a 
sequel yet. 

Particle tracking experiments in a colloidal suspension of PMMA particles 
revealed an effective temperature of the order of double the ambient one from 
the mobility-diffusivity relation [93] . 

In the soft matter realm a favorite is an aqueous suspension of clay, Laponite 
RG, in its colloidal glass phase. During aging, because of electrostatic attraction 
and repulsion, Laponite particles form a house-of-cards-like structure. After 
a number of rather confusing reports the status of Teff in this system can 
be summarized as follows. The surprisingly high Teff found with dielectric 
spectroscopy combined with spontaneous polarization noise measurements was 
later ascribed to violent and intermittent events possibly linked to the presence 
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of ions in the solution which may be the actual source of FDT violation. For 
the moment dielectric degrees of freedom are invalidated as a good test ground 
for Teff in this sample [94] . Using other methods several groups found that Teff 
detaches from the bath temperature. Strachan et al. [95] measured the diffusion 
of immersed probe particles of different sizes via dynamic light scattering and 
simultaneous rheological experiments and found a slightly higher Teff than T . 
With micro-rheology Abou and Gallet observed that Teff increases in time from 
T to a maximum and then decreases back to T [96]. Using a passive micro- 
rheology technique and extracting Teff from the energy of the probe particle 
via equipartition Greinert et al. also observed that Teff increases in time [97]. 
In parallel, a series of global mechanical tests, and passive and active micro- 
rhcological measurements that monitor the displacement and mobility of probe 
Brownian particles were performed by Ciliberto's and D. Bonn's groups, both 
finding no violation of FDT over a relatively wide frequency range [94, 98, 99]. 
In a very detailed article Jop et al. explain many subtleties in the experimental 
techniques employed and, especially, the data analysis used to extract Teff that 
could have biased the results quoted above. A plausible reason for the lack of 
out of equilibrium signal in some experiments using Laponite as well as other 
colloidal glasses is that the range of frequency-time explored may not enter the 
aging regime. Moreover, none of these works studied the degrees of freedom 
of the Laponite disks themselves but, instead, the properties of the solvent 
molecules or probe particles. More recently, Maggi et al. combined dynamic 
light scattering measurements of the correlation function of the colloid rotations 
with those of the rcfringence response [100] and a x(C, ^tu) plot that is rather 
constant as a function of C and slowly recovers the equilibrium form as the 
arrested phase is approached (t^, ranges from 90 to 1200 min and the violations 
arc observed for time differences between 0.1 and 1 ms, i.e. frequencies between 
10 and 1 kHz). Teff is at most a factor of 5 larger than T. The actual behaviour 
of Laponite remains mysterious - and not only in what Teff is concerned! 

Oukris and Israeloff measured local dielectric response and polarization noise 
in polyvinyl-acetate with electric-force-microscopy [101]. They probed long- 
lived nano-scale fluctuations just below Tg, achieved a good signal-to- noise ratio 
down to very low frequencies, constructed a parametric plot by keeping t^ fixed 
and found a non-trivial asymptotic form with no t^ dependence within the avail- 
able accuracy. The data combine into the parametric plot Teff(C') ~ 2^(7-0.57 
in the aging regime. 

5.6 Relaxation in frustrated magnetic systems 

Disordered and frustrated magnets behave collectively at low temperatures and 
developed ordered phases that although not fully understood are accepted to 
exist. As macroscopic glassy systems they present a separation of time-scales in 
their low-temperature dynamics and are good candidates to admit a thermody- 
namic interpretation of the FDT violations. 
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5.6.1 Remarks on model systems 

The physics of spin-glasses is a controversial subject. Some authors push an Ising 
domain-growth interpretation of their dynamics - slowed down by domain wall 
pinning by disorder - a.k.a. the droplet picture [102]. If the scheme discussed 
in Sect. 5.2 were reproduced under strong disorder, the asymptotic x(C') plot 
would have a linear piece of slope — l/T and a sharp transition at to a flat 
aging piece. The domain-growth interpretation is not accepted by other authors 
and more complex scenarii based on the static [30] and dynamic [6] solution to 
the SK model are envisaged, with a non-trivial x(C) as a result. Much effort has 
been put in trying to interpret numerical and experimental data as validating 
one description at the expense of the other. Unfortunately, it is very difficult 
to distinguish between the two. A third possibility is that, in a loose sense, the 
spin-glass be like the low-T phase in the 2d XY model, with quasi long-range 
order. Yet another proposal is that actual spin-glass samples are of Heisenberg- 
type and that chilarity might be decoupled from spin with a chiral-glass order 
arriving at a higher critical temperature than the spin-glass ordering [103]. 

The trap model [104] was devised to describe slow dynamics in systems with 
weak crgodicity breaking and it was applied, notably, to describe experiments 
in spin-glasses. The model shows a glass transition at a Tg below which an 
equilibrium Boltzmann state cannot exist. The xi^) ^ slope that varies 
continuously even though there is a single scaling of relaxation times with age, 
it depends non-trivially on the observable and one cannot use it to define a 
meaningful Teff [38]. The reason for this failure seems to be the unbounded 
nature of the energy and the fact that an equilibrium distribution does not 
exist below Tg. 

5.6.2 Simulations 

Monte Carlo simulations of the 3d Edwards- Anderson (EA) model were carried 
out by several groups. One of the hallmarks of the dynamics of the SK model, 
dynamic ultrametricity [6] , is absent from all numerical and experimental data 
analyzed so far. Magnetic correlation and susceptibility relax in two scales, 
the by now usual stationary one for finite time-differences and an aging one in 
which the data are well described by a simple t/tw scaling. This aging scaling 
does not conform with the droplet picture either, which predicts an asymptotic 
Int/lntw form. In all studies so far, the parametric plot was constructed by 
keeping t^ fixed and the curves drift towards increasing values of x for longer 
tu,s as in a transient or critical system. In simple coarsening problems the drift 
with increasing t^ goes in the opposite direction of rendering the aging part of 
the curves flatter; this remark suggests to discard a simple droplet picture. The 
outcome x(C) found for the longest t^^ reached was interpreted as being non- 
constant [105] - as in the SK model - although this is, in our opinion, not that 
clear from the data that could be described by a straight line. The simultaneous 
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t/tw scaling, the lack of unambiguous evidence for a stable plateau at Qea, and a 
curved x(C) in the aging regime is not what would be expected from an analogy 
with the SK model. Instead, it would be consistent with critical dynamics and 
the 2d XY model similitude. A number of caveats on the numerical analysis 
should, however, be lifted before reaching a firm conclusion. 

The finite-time finite-length relation between static H(C, ^(itu)) and long- 
time out of equilibrium dynamic x(C, ^u,) (Sect. 3.5.1) [41] was put to the test 
in the 2d and ?>d EA models at finite T. The notable coincidence of the two 
functions found in the 2d case, in which there is no complex equilibrium struc- 
ture, suggests that the claimed coincidence of x(C) and H(C) in Sd [105] might 
also be valid just in the transient regime. 

Simulations of the id Heisenberg spin-glass model with weak anisotropy 
suggest that Teff associated to the spin degrees of freedom is constant and 
about twice the critical temperature for spin-glass ordering [106]. As far as we 
know, chiral degrees of freedom have not been used to estimate Teff- 

As regards fluctuations, the two kinds discussed in Sect. 4.7 were measured 
in the 3d EA spin-glass. Disordered induced ones [52], in which one computes 
strictly local noise-averaged correlations and linear responses, demonstrate the 
existence of two types of spins in each sample: rapid paramagnetic-like ones and 
slow ones. The former satisfy FDT while the latter evolve in two time-regimes 
with a fast one satisfying FDT and a slow one in which Xi(Ci) looks quite flat as 
in coarsening systems. The simulation suggests that the two ensembles behave 
independently of each other and are strongly correlated with the backbone of 
the ground state configurations. The average over all sites (at finite t^) gives 
rise to a curve with non-constant slope. These results suggest a still different 
picture for the spin-glass dynamics in which a rather compact set of spins un- 
dergoes coarsening of the backbone equilibrium configurations while the other 
ones behave paramagnctically. This intriguing idea needs to be put to further 
test. 

The analysis of noise induced fluctuations suggests that cq. (36) is valid 
although better numerical data would be needed to have definitive evidence 
for this statement. A more detailed discussion can be found in [56]. Very 
recent studies of non-linear fiuctuations that take advantage of FDRs to compute 
higher order responses point in the direction of the TRI scenario [59] with a finite 



5.6.3 Experiments 

On the experimental side the first attempt to quantify FDT violations in spin- 
glasses was indirect [107]. Simultaneous measurements of global magnetic noise 
and susceptibility in the thiospinel insulating spin-glass were later performed 
by Herisson and Ocio [108]. The data confirm deviations from the FDT with 
a x(C, *■!«) plot of relatively curved form although still evolving during the ex- 
perimental time window. The authors interpreted it as evidence for the full 
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RSB scenario, via the association x(C) ^ ^(C*)- However, as with numerical 
data [105], dynamic ultrametricity fails to show off, the asymptotic limit of the 
parametric construction is still far, and a clear-cut distinction between a curved 
and a linear x(C) is hard to assess. 

More recent experiments exploit two novel techniques, Hall-scnsor based 
magnetometer and giant magnetoresistance technology to detect signals from 
very small samples [109]. The use of these probes opens the way to perform 
a systematic study of FDT violations in magnetic systems of different kind 
(spin-glasses, super-spin glasses, disordered ferromagnets...). The first of these 
measurements appeared recently [110] in a super-spin glass, a system of mag- 
netic nanoparticles suspended in fluid glycerol with a single-domain magnetic 
structure that behaves as one large spin, the orientation of which is the only 
degree of freedom. The large magnetic moment facilitates the observation of 
magnetic noise. For aging times of the order of 1 h, the ratio of Teff to the 
bath temperature T grows from 1 to 6.5 when T is lowered from Tg to 0.3 Tg, 
regardless of the noise frequency. 

Artificial spin ice is yet another material in which the Teff notion has been 
tested [111]. 

5.7 Driven liquids and glasses 

In [112] the molecular dynamics of a binary Lennard- Jones mixture under a 
steady and homogeneous shear flow was studied. The deviation from FDT 
is similar to the one found analytically in disordered spin models of RFOT 
type with asymmetric couplings that mimic non-conservative forces [57, 113]. 
Moreover, it does not depend on the observable. The tracer particle experiment 
was also realized. When the tracers' Einstein frequency is smaller than the 
inverse relaxation time of the fluid, a non-equilibrium equipartition theorem 
holds with mTRW? = Teff, where Vz is the velocity in the direction transverse 
to the flow. For increasing the effective temperature very slowly crosses 
over from T to the slow modes value, in perfect agreement with the notion 
of a temperature measured by a thermometer sensible to the scale. Teff also 
captures the essential phenomenological idea that when a system is sheared 
more vigorously its effective temperature increases. 

O'Hern et al. also studied fluctuation-dissipation relations in shear flu- 
ids [76]. This group defined an effective temperature through the 'static limit' 
limf^oo x{t ~ iw)/C{t,t), a kind of average of the slope of the x(C) plot over 
the full range of C(i — t^) that mixes different time scales (in particular, the 
high and low frequency ones) . A more thorough discussion of the comparison 
between this definition and the one described in this review was given by Ilg 
and Barrat [75] within a fully solvable model that demonstrates the importance 
of not mixing time-scales to get physically sensible results. 

A first study of the fiuctuations of entropy production in a Lennard- Jones 
fluid above and below Tg under a shear fiow appeared in [114] and the need to 
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take into account Teff, as obtained from the modification of the FDT below Tg, 
was signaled in this paper. A more detailed analysis of the time-scale dependent 
effective temperature would be needed to fully test the proposal in [49] . 

Another prominent example is the current driven motion of vortices in 
type II superconductors. Disorder reduces dissipation, is responsible for non- 
equilibrium transport and magnetic properties. The external force induces two 
dynamic phase transitions separating plastic flow, smectic flow and a frozen 
transverse solid. A low-frequencies Teff that decreases with increasing driv- 
ing force and reaches the equilibrium melting temperature when the dynamic 
transverse freezing occurs was computed from the transverse motion in the fluid 
moving phase [115]. 

5.8 Granular matter 

Several studies of the effective temperature of granular matter have been pur- 
sued theoretically, numerically and experimentally. In the latter front, D'Anna 
et al. immersed a torsion oscillator in a granular system fluidized by strong high 
frequency external vibrations to realize the 'thermometer' experiment. They 
found Teff cx with r the adimensional measure of vibrational intensity, and 
quite independently of w [116]. Wang et al. [117] visualized the dynamics of 
tracer particles embedded in a id granular ensemble slowly sheared by the ro- 
tating inner wall of a Couette cell. Teff, as obtained from the comparison 
between the tracer's diffusion and mobility perpendicular to the applied rate 
of strain, is independent of the shear rate used and the tracers properties but 
does depend on the packing density of the system. Tests of the thermodynamic 
properties of Teff have not been carried through in this system yet. The depen- 
dence on the direction of the applied stress was studied by Twardos and Dennin 
in a plastic bead raft close to jamming [118]. As expected, the correlations 
and linear responses in the direction of flow do not decay slowly and does 
not have the same properties as in the transverse direction (cfr. [119] and [69]). 
Gei and Behringer stressed the fact that in a granular assembly the outcome 
of a mobility measurement depends on whether one imposes the velocity or the 
external force [120]. 

In the powders literature reference is often made to the 'granular tempera- 
ture', a measure of the temperature of the fast modes, as given by the kinetic 
energy of the grains = ji^K = Importantly enough, Tr- is a high 

frequency measure that does not really access the structural properties of the 
sample and, in a sense, plays the role of the environmental temperature in ther- 
mal systems. Tk is generically smaller than Teff, as in thermal systems where 
Tk = T, the temperature of the bath. 
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5.9 Activated dynamics 

Activated processes often occur in systems that are out of equilibrium, in the 
sense that their response to an external drive is strongly non-linear or that 
their phase space distribution is not the Gibbs-Boltzmann one. The question 
as to whether an Arrhenius law governs the activation rate, possibly with an 
effective temperature, and how the latter compares to the one defined from the 
deviations from FDT has been addressed recently [121]. Ilg and Barrat studied 
the effect of an out of equilibrium flowing environment, a weakly sheared super- 
cooled liquid, on the activated dynamics between the two stable conformations 
of dumbbell particles. The transition rate is well described by an Arrhenius 
law with a temperature that crosses over from the one of the equilibrium bath 
to a higher value close to the Teff of the slow modes of the driven fluid. The 
crossover roughly occurs at the value of the rate that corresponds to the inverse 
of the a relaxation time of the fluid. 

Three related studies are also worth mentioning. An effective temperature, 
also consistent with the one stemming from fluctuation-dissipation measure- 
ments, appears in a phenomenological Arrhenius law that describes transverse 
jumps between channels in the driven motion of vortex lattices with random 
pinning [119]. Haxton and Liu showed that in the shear dominated regime the 
stress of a 2d sheared fluid follows an Arrhenius law with the effective temper- 
ature [122]. A study of activation and Teff in a 2d granular system close to 
jamming appeared in [123]. 

5.10 Biological systems 

In biologically inspired problems the relevance of Teff was stressed to reveal the 
active process in hair bundles [125] and model cells [126]. Morozov et al. [127] 
studied a model of the cytoskeletal network made of semi-flexible polymers sub- 
ject to thermal and motor-induced fluctuations and found a Teff that exceeds 
the environmental temperature T only in the low-frequency domain where mo- 
tor agitation prevails over thermal fluctuations. Simple gene network models 
were studied from the Teff perspective in [124]. Fluctuation-dissipation ratios 
were used to quantify the degree of frustration, due to the existence of many 
metastable disordered states, in the formation of viral capsids and the crys- 
tallization of sticky discs, two self-assembly processes [128]. Fluctuations and 
responses of blood cell membranes for varying ATP concentration were mea- 
sured very recently [129]. The measured Teff approaches the bath temperature 
at high frequencies and increases at low frequencies reaching 4-10 times the 
ambient temperature. 

Ratchets are simple models of molecular motors, out of equilibrium systems 
with directed dissipative transport in the absence of any external bias. Harada 
and Sasa proposed to use the violations of FDT in flashing ratchets as a means 
to measure the energy input per unit time in molecular motors - an otherwise 
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difficult quantity to access [130]. Kolton showed that the rectified transverse 
velocity of a driven particle in a geometric ratchet is equivalent to the response of 
a Id flashing ratchet at a drive-dependent Teff, as defined from the generalized 
Einstein relation [131]. 

Active matter is driven out of equilibrium by internal or external energy 
sources. Its constituents absorb energy from their environment or from inter- 
nal fuel tanks and dissipate it by carrying out internal movements that lead 
to translational or rotational motion. A typical example are self-propelled par- 
ticle assemblies in bacterial colonies. The role played by Teff in the stabil- 
ity of dynamic phases of motorized particle systems was stressed by Shen and 
Wolynes [132]. Multiple measurements of Teff were carried out with molecular 
dynamic simulations of motorized spherical as well as linear molecules in in- 
teraction [133]. All measurements (from fiuctuation-dissipation ratio and using 
tracers) yield a constant low-frequency Teff > T when the effect of the motors 
is not correlated with the structural rearrangements they induce. Instead, Teff 
takes a slightly lower value than T when susceptible motors are used, as argued 
in [132]. Such an 'inversion' also occurs in relaxational systems in which the 
initial configuration is chosen to be one of equilibrium at a lower T than the 
working one [43, 44, 45]. In the case of uncorrelated motors, Teff/T was found 
to follow the empirical law Teff/T ~ 1 + 7/^ with / the active force relative to 
the mean potential force and 7 ^ 15 a parameter. Palacci et al. [134] inves- 
tigated Teff by following Perrin's analysis of the density profile in the steady 
state of an active colloidal suspension under gravity. The active particles used 
- JANUS particles - are chemically powered colloids and the suspension was 
studied with optical microscopy. The measurements show that the active col- 
loids arc hotter than in the passive limit with a Teff that increases as the square 
of the parameter that controls activation, the Pcclct number, a dependence that 
is highly reminiscent of the /^-dependence of the simulations mentioned above. 
Other theoretical studies of Teff in active matter appeared in [135]. 

Joly et al. [136] used numerical techniques to study the non-equilibrium 
steady state dynamics of a heated crystalline nanoparticle suspended in a fluid. 
This problem models an active colloid that acts as a local heat source and 
generates a temperature gradient around it. By comparing the mobility to the 
velocity correlation function, they found that the FDT approximately holds at 
short-time lags with a temperature value that coincides with the kinetic one. In 
contrast, at long-time lags data arc compatible with the temperature estimated 
by using the Einstein relation. 

Certainly, many more studies of effective temperatures will appear in this 
very active field of research, essentially out of equilibrium, in the near future. 

5.11 Plasticity 

Langer et al. extended the traditional phenomenological defect-fiow theory to 
the shear transformation zone (STZ) theory of large-scale plastic deformation 
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in amorphous materials [137]. The new theory incorporates effective tempera- 
ture ideas. It is a picture of plastic deformation in molecular glasses in which a 
'disorder temperature' characterizes the steady state of the system and controls 
the slow processes. The relation between the disorder temperature and a con- 
figurational entropy, under the assumption of a sharp separation of time scales 
between structural and vibrational processes, as well as other thermodynamic 
properties of it - along the lines discussed in Sect. 4.8 - were addressed. As 
far as we know, there have been no tests to compare the STZ disorder tem- 
perature and the effective temperature computed from fluctuation-dissipation 
measurements; this is the reason why we keep distinct names for the two quan- 
tities. The STZ theory suggests an explanation of shear-banding instabilities 
that have been put to the numerical test in [138]. 

Another approach to plasticity consist in adapting the RFOT-replica ap- 
proach to this problem [139]. Numerical simulations of binary soft sphere mix- 
tures at low temperature prove that the stress relaxation and its response to a 
strain step are also linked by a modified FDT with a single valued Teff in the 
aging regime [140] . 

5.12 Turbulent fluids 

Last but not least, the experimental quest for effective temperatures in the 
steady state of turbulent flows has recently restarted. As far as we know, two 
experiments on turbulent flows appeared in the literature after the role played by 
different time-scales was stressed in the investigation of effective temperatures 
in macroscopic systems with slow dynamics. 

The transverse fluctuations and the response of a string held at its ends and 
at constant tension (a mechanical probe) in the inertial regime of a station- 
ary turbulent air jet flow were compared to obtain the effective temperature. 
The ranges of wave-numbers, and Reynolds number (7.4 10^ < Re < 1.7 10^). 
accessed in this experiment are pretty wide [141]. All measurements are compat- 
ible with Teff oc fc~^^/^. The k dependence conflrms that there is no equilibrium 
between Fourier modes due to the energy flux between scales. The particular 
exponent is explained in [141] with a simple model derived from Kolmogorov 
1941. 

In a more recent experiment focus was set on an anisotropic non-homogeneous 
axisymmctric von Karman flow at large Reynolds number [142]. The measure- 
ments are, though, somehow indirect and yield an effective temperature that 
depends on the observable. This feature needs to be better characterized. 

5.13 Quantum models 

There is growing interest in the dynamics of quantum systems. Of particular 
importance in this field is to distinguish cases in which the system of interest 
is isolated and the dynamics occur at constant energy from those in which the 
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system is coupled to an environment and the dynamics are dissipative. The out 
of equihbrium dynamics are typically induced in two ways: the system is driven 
out of equilibrium by, for instance, a coupling to electric or heat current sources, 
or it is strongly perturbed by time-dependent external fields; the system evolves 
after a quench meaning that a parameter in its Hamiltonian is changed with 
some protocol. All these cases are easily realized in the laboratory nowadays 
and their potential applicability is being explored. Theoretical searches for 
effective temperatures in these problems, defined in different ways, are starting 
to appear in the literature. 

As recalled in Sect. 2 linear responses and correlation functions in equilib- 
rium quantum systems are related by FDTs that take slightly different form 
depending on the bosonic or fermionic character of the observables. In an out 
of equilibrium situation one can compare the linear response and the correlation 
and test whether, at least in some dynamic regime, a parameter replacing the 
environmental temperature appears. This route was first taken within mean- 
field quantum glassy models quenched from their disordered into their ordering 
phase [36, 37, 143]. In these models, the quantum FDT holds within the rapid 
stationary scale. Instead, when the dynamics gets slow and the aging regime is 
attained the relation between linear response and correlation takes the classical 
form, with an effective temperature that is higher than the one of the bath and 
different from zero. As in the classical case, the static properties of these models 
can be solved with the replica and TAP approaches and a connection with the 
dynamic solution can be established [144] . 

More recently, the analysis of mesoscopic quantum models commenced. A 
metallic ring threaded by a time-dependent magnetic field and coupled to a 
lead in equilibrium at temperature T and chemical potential fi was studied 
in [145]. The numerical solution to the Schwingcr-Keldysh equations for the 
free but driven fermions, in the limit of small dissipation, suggests that their 
Green functions satisfy an FDT with constant Teff. More recently, Arrachea 
et al. attacked the problem of a wire connected to left and right reservoirs 
(in equilibrium at the same chemical potential and temperature) and driven 
out of equilibrium by different ac pumps locally connected to the wire [146]. 
The local effective temperature was computed from the modified FDT and by 
requiring that there be no heat flow to nor from macroscopic probes, i.e. ther- 
mometers, weakly coupled to chosen sites on the device. For weak driving and 
environmental temperatures lower than the Fermi energy of the electrons these 
two measurements coincide on each site. Teff has spatial 2fcj;^-Friedel-like os- 
cillations (allowing for local cooling) but its spatial average is higher than the 
temperature of the reservoirs when the pump is applied. Moreover, the direc- 
tion of heat flow between the device and each one of the leads is dictated by 
the temperature gradient at the contact, defined as the difference between Teff 
at the contact and the temperature of the lead. For conveniently chosen pumps 
the device can extract energy from one lead and transport it to the opposite 
one. 
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The theory of isolated quantum models is a very active field of research, 
pushed by the large number of experiments on cold-atom systems that are being 
realized World-wide. Questions on thermalization after a quantum quench, in 
which a parameter in the Hamiltonian is suddenly changed, are being posed, 
extending this very much investigated issue in classical statistical mechanics 
to quantum systems. At present the belief is that non-integrable quantum 
systems do reach equilibrium while integrable ones do not. Rigol et al. proposed 
to use extended Gibbs-Boltzmann measures in the latter case with as many 
constants of motion as necessary fixed by the initial condition [147]. Currently, 
this proposal is been debated and studied in particular cases. For instance, 
Rossini et al. complemented it with the idea that non-local observables (when 
expressed in terms of quasi-particles) should equilibrate in integrable systems 
while local ones should not under general conditions [148]. They based this 
conjecture on their mixed analytic/numeric solution of a quantum spin chain. 
Studies of effective temperatures (as defined from deviations from FDT) in this 
context arc under way [149]. Another quantum setting in which an effective 
temperature was proposed is the entanglement entropy of a sub-system [150]. 

6 Conclusions 

In this review we discussed the notion of an effective temperature [l]-[4] that 
stemmed from the deviations from the fiuctuation-dissipation theorem found in 
the analytic solutions to simple glass models [5, 6] and self-consistent approx- 
imations to more realistic glassy systems [32, 33]. Although phenomenological 
definitions of out of equilibrium temperatures were not lacking in the literature, 
see e.g. [8, 10], the effective temperature discussed in this report lies, in our 
opinion, on a firmer physical basis. 

Fluctuation dissipation violations can be searched for in all out of equilibrium 
system. Whether the outcome can be interpreted as giving origin to an effective 
temperature is a different and more delicate issue. In this review we tried to 
distinguish cases in which the latter can be done from cases in which it is not 
possible. A necessary condition seems to be that the system should evolve 
slowly, irrespectively of whether it relaxes or is in a non-equilibrium steady 
state, in a small entropy production regime [4]. Although not rigorously proven 
it seems natural to require that it approaches an approximately flat measure on 
a region of phase space, the entropy of which (also called complexity) should 
give an alternative, microcanonic-like, access to the effective temperature [2, 7, 
21, 63]. A consequence of the last condition is that the energy density, and 
more generally the averages of one-time dependent observables, converge to 
finite values. All these features are satisfied exactly in some solvable mean- 
field-like glasses such as the celebrated p-spin disordered model and the low-T 
mode-coupling approximations. Numerical simulations suggest they also are, at 
least within a given time regime, in a large number of glassy systems with short- 
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range interactions, including Lennard- Jones mixtures and others. In all these 
cases the systems have a rather complex collective behavior with a separation 
of time-scales fast-slow that has an influence in the values that the effective 
temperature takes. 

As far as it has been checked, in all the above-mentioned cases the effec- 
tive temperature conforms to the common prejudices one has of a temperature 
and, more importantly, it has the most welcome property of being measurable 
directly, hence being open to straightforward ~ though difficult to implement 
- experimental tests. Central to the correct identification is the realization 
that the relaxation time-scales have to be correctly identified and that measure- 
ments have to focus on each of them separately. Simulations in atomic glasses 
are very complete and yield pretty convincing support to the effective temper- 
ature ideas [85]-[89],[112]-[114]. Experimental subtleties prove to be difficult to 
keep under control, especially in colloidal suspensions [93]- [100], but successful 
measurements of fluctuation-dissipation violations in structural glasses [92, 101] 
have been surveyed. Very recent numerical [133] and experimental [134] studies 
of active matter are also very promising. Simulations of spin-glass [105] clearly 
demonstrate the violation of FDT but are not as precise as to determine be- 
yond doubt the actual functional form of the modified relation. In our opinion, 
pre-asymptotic effects are still to be disentangled from the asymptotic regime 
in a satisfactory way. The same applies to the first experiments in an insulator 
spin-glass [108] the outcome of which could be greatly improved by the use of 
innovative techniques [109, 110]. A careful analysis of the thermodynamic prop- 
erties of the fluctuation-dissipation ratio remains to be done such frustrated and 
disordered magnets. 

Quite naturally, the idea that effective temperatures could be relevant to 
other out of equilibrium systems was explored in the last flfteen years. Crit- 
ical quenches in dissipativc classical systems [12], zero temperature dynamics 
at the lower critical dimension [13], granular matter [68]-[73] and turbulent flu- 
ids [141, 142] are just some examples in which fluctuation dissipation relations 
have been studied. The question remains as to whether the outcome admits a 
thermodynamic interpretation in, at least, some limits. In the case of critical 
quenches, this seems to be the case when the long two-time limit prescription 
of Godreche and Luck is taken [79]. Simulations of weakly perturbed granu- 
lar matter are very encouraging [69] but experiments have turned out to be 
much harder to realize [67] (see, though, [117]). In turbulent fluids the question 
remains open experimentally. 

Out of equilibrium quantum systems are receiving enormous attention nowa- 
days, boosted by experiments in cold atoms and nanotechnology. In this realm, 
fundamental questions of thcrmalization arise, especially in isolated samples 
submitted to a quench. The effective temperature has been studied in a few 
dissipative quantum systems out of equilibrium, both mean-field [36, 37] and 
low-dimensional [145, 146]. The analysis of whether it also plays a role in quan- 
tum quenches of isolated samples is under way [149]. 
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This article presents a vast panorama of fluctuation dissipation deviations in 
non-equilibrium classical, and to a much smaller extent, quantum systems, and 
their interpretation, in some cases, in terms of effective temperatures. Many 
questions of fundamental interest remain open in this field. One of the main 
open challenges in the context of glassy systems, is to find the microscopic origin 
of these modifications. From a wider viewpoint the validity of the effective 
temperature concept should be more clearly delimited. 
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